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Applications of perturbative QCD to deeply virtual Compton scattering and hard exclusive elec- 
troproduction processes require a generalization of usual parton distributions for the case when 
long-distance information is accumulated in nonforward matrix elements {p'\O(0, z)\p) of quark 
and gluon light-cone operators. We describe two types of nonperturbative functions parametriz- 
ing such matrix elements: double distributions F(x, y; t) and nonforward distribution functions 
!F((X;t), discuss their spectral properties, evolution equations which they satisfy, basic uses and 
general aspects of factorization for hard exclusive processes. 
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I. INTRODUCTION 



The standard feature of applications of perturbative QCD to hard processes is the introduction of phenomenolog- 
ical functions accumulating information about nonperturbative long-distance dynamics. The well-known examples 
are the parton distribution functions f p /H {%) used in perturbative QCD approaches to hard inclusive processes 
and distribution amplitudes ip^(x), lpn(xi, %2, £3), which naturally emerge in the asymptotic QCD analyses of 
hard exclusive processes [|[j7). Recently it was argued that the same gluon distribution function f g (x) used for 
description of hard inclusive processes also determines the amplitudes of hard exclusive J/ip || and p-meson |9| 
electroproduction. Furthermore, it was proposed |]lof to use another exclusive process of deeply virtual Compton 
scattering 7*p — > 77/ (DVCS) for measuring quark distribution functions inaccessible in inclusive measurements 
(earlier discussions of nonforward Compton-like amplitudes ^*p — > j*p' with a virtual photon or Z° in the final 
state can be found in refs. p|-p^|). The important feature (noticed long ago ]l]]Jl^| ) is that kinematics of hard 
elastic electroproduction processes (DVCS can be also treated as one of them) requires the presence of the longi- 
tudinal component in the momentum transfer r = p — p' from the initial hadron to the final: rn = CP- F° r DVCS 
and p-electroproduction in the region Q 2 » \t\,m 2 H , the longitudinal momentum asymmetry (or "skewedness" ) 
parameter Q coincides with the Bjorken variable XBj = Q 2 /2(pq) associated with the virtual photon momentum q 
[pr[ . This means that the nonperturbative matrix element (p' \ . . .\p) is essentially asymmetric and, strictly speak- 
ing, the distributions which appear in the hard elastic electroproduction amplitudes differ from those studied in 
inclusive processes. In the latter case, one has a symmetric situation when the same momentum p appears in both 
brackets of the hadronic matrix element (p\ . . . \p). 

For diffractive processes, one deals with a kinematic situation when both the variable £ specifying the longitudi- 
nal momentum asymmetry (skewedness) of the nonperturbative matrix element (p'\ . . . \p) and the absolute value 
of the momentum transfer t = (p 1 — p) 2 are small. Studying the DVCS process, one should be able to consider the 
whole region < £ < 1 and t ~ 1 GeV 2 In this situation, one deals with essentially nonforward (or off-forward 
in terminology of Ref. jio]]) kinematics for the matrix element (p'\ . . . \p). The basics of the PQCD approaches incor- 
porating asymmetric/off-forward parton distributions were formulated in refs. |l0| , pl|p^ , |i"5[ . A detailed analysis of 
PQCD factorization for hard meson electroproduction processes was given in Ref. f [f8[. Applications of asymmetric 
gluon distributions to elastic diffractive J/ip electroproduction were discussed in ]l9|-^l| . In a recent paper (22) , 
the off-forward quark distributions were studied within the MIT bag model. A numerical study of the evolution of 
the asymmetric gluon distribution was attempted in Ref. |p3[ . Thus, there is an increasing interest in the studies 
of these new types of hadron distributions, their general properties and applications. 

Our goal in the present paper is to give a detailed description of the approach outlined in our earlier papers 



17 1 . The basic idea of refs. ]1^,0 is that constructing a consistent PQCD approach for amplitudes of hard 
exclusive electroproduction processes one should treat the initial momentum p and the longitudinal part of the 
momentum transfer r on equal footing by introducing double distributions F(x 1 y), which specify the fractions 
of p and r- 1| carried by the constituents of the nucleon. These distributions have hybrid properties: they look 
like distribution functions with respect to x and like distribution amplitudes with respect to y. Writing matrix 
elements of composite operators in terms of double distributions is the starting point of constructing the PQCD 
parton picture. Another important step is taking into account the logarithmic scaling violation. The evolution 
kernels R(x,y; £,77) for double distributions have a remarkable property: they produce the DGLAP evolution 
kernels P(x/£f) ]2^-|2^| when integrated over y, while integrating R(x, y; £, rf) over x one obtains the BL-type 
evolution kernels V{y, rf) HQ for the relevant distribution amplitudes^ One can use these properties of the kernels 
to construct formal solutions of the one- loop evolution equations for the double distributions. The longitudinal 
momentum transfer r\\ is proportional to p: rn = C^p and, for this reason, it is convenient to parametrize matrix 
elements (p — r \ . . . \p) by asymmetric distribution functions Tq(X) specifying the total light-cone fractions Xp, 



^Originally, the evolution equation for the pion distribution amplitude in QCD was derived and solved in Ref. ^), where 
the anomalous dimension matrix Z n k was used instead of V(y,n) (see also |27|). 
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(X — £)p of the initial hadron momentum p carried by the "outgoing" and "returning" partons^. It should be 
emphasized that double distributions F(x, y) are universal functions in the sense that they do not depend on the 
skewedness parameter £ while the asymmetric distribution functions T^{X) form a family of X-dcpcndcnt functions 
changing their shape when ( is changed. The functions !F^(X) also have hybrid properties. In the region X > ( 
the returning parton carries a positive fraction (X — £)p of the initial momentum p, and hence J-^(X) is similar to 
the usual parton distribution f(X). On the other hand, in the region < X < £ the difference X — £ is negative, 
i.e., the second parton should be treated as propagating together with the first one. The partons in this case share 
the longitudinal momentum transfer r|| = Qp in fractions Y = X/C, and 1 — Y. This means that in the region 
X < ( the function T^{X) looks like a distribution amplitude. It is possible to formulate equations governing the 
evolution of the asymmetric distribution functions !Ff(X) and establish relations between these functions, double 
distributions F(x,y) and usual distribution functions f(x) 16 17[] . 



Constructing a QCD parton picture for hard electroproduction processes, it is very important to know spectral 
properties of the relevant parton distributions F(x,y) and ^(X). Using the approach |2^] based on the a- 
representation analysis [p9|"p2[, it is possible to prove that double distributions F(x,y) have a natural property 
that both x and y satisfy the "parton" constraints 0<jc<1,0<j/<1 for any Feynman diagram contributing 
to F(x,y). A less obvious restriction < x + y < 1 guarantees that the argument X = x + yC of the asymmetric 
distribution J-^(X) also changes between the limits and 1. An important observation here is that X = can 
be obtained only if both x — and y — 0. Because of vanishing phase space for such a configuration, one may 
expect that asymmetric distributions T^X) vanish for X — 0. This property is very essential, because the hard 
subprocess amplitudes usually contain l/X factors. When F^(0) ^ 0, one faces a singularity Fq(X)/X at the 
end-point of the integration region < X < 1. Since such a singularity is not integrable, factorization of short- 
and long-distance contributions does not work in that case. 

The paper is organized in the following way. In Sec. II, we consider parton distributions in a toy scalar model. 
Despite its simplicity, it shares many common features with the realistic QCD case. In particular, the spectral 
properties of distribution functions are not affected by the numerators of quark and gluon propagators, derivatives 
in triple-gluon vertices, etc. Hence, studying a scalar model we just concentrate on the denominator structure of 
the relevant momentum integrals, which is the same in both theories. We start with the simplest example of the 
usual (forward) distribution f(x) and then consider more and more complicated functions: the double distribution 
F(x,y), asymmetric distribution function Tq(X) and nonforward distribution !Ft(X\t). Explicit expressions for 
these functions at one-loop level are obtained with the help of the a-representation. Using the latter one can 
easily establish the spectral properties of the distribution functions. The a-representation also provides a very 
effective starting point for a general analysis of factorization and large-Q 2 behavior of elastic amplitudes. In 
Sec. Ill, we outline the all-order extension of the one-loop analysis. We give an all-order definition of the double 
distribution function F{x,y) and demonstrate that it has the spectral properties < {x, y,x + y} < 1. We show 
how one can use the a-reprcscntation analysis for finding integration regions responsible for the leading large-Q 2 
contributions. We also discuss modifications of twist counting rules in QCD due to cancellations between different 
gluonic contributions in the Feynman gauge and other complications which appear in gauge theories. In Sec. IV, 
we give definitions of nonforward distributions Tq{X] t) in QCD. Just like the usual distribution functions f{x) and 
distribution amplitudes <p(y) 7 the new distributions depend on the factorization scale \i, i.e., it is more appropriate 
to use the notation Tc(X\ t; p) for the nonforward distributions rather than simply T(^{X\ t). Evolution equations 
governing the /x-dependence of the nonforward distributions are discussed in Sec. V and Sec. VI. We show how 
one can obtain evolution kernels for nonforward distributions using already known kernels B{u 1 v) of the evolution 
equation for the light-ray operators J33|,[54| . Since this equation has an operator form, substituting it into a specific 
matrix element one can convert B{u, v) into desired evolution kernels. In particular, taking (p\ . . . \p) one obtains 
the DGLAP kernels, choosing (0| . . . \p) one gets BL-type kernels while resorting to (p' \ . . . \p) and parametrizing the 
matrix elements through F(x,y) or Tq(X) one ends up with the kernels R(x,y; £,rj) and W^(X,Z) governing the 



*The asymmetric distribution functions defined in Ref. |17j are similar to, but not coinciding with the t — ► limit of the 
off- forward parton distributions introduced by X.Ji [hoi, see Sec. IX. 



evolution of double and asymmetric/nonforward distributions, respectively. In Sec. V, we discuss the derivation 
of the evolution kernels Wq(X,Z) for the nonforward distributions. We show, in particular, that in the region 
< {X,Z} < C, the kernels W^(X,Z) reduce to the BL-type kernels V(X,Z) calculated for rescaled variables 
X/(, Z/Q. This result is very natural, since !F^{X) can be treated as a distribution amplitude when X < £. In 
the opposite limit ( < {X, Z} < 1, the evolution is similar to that of the DGLAP equation, the basic distinction 
being the difference between the outgoing X, Z and returning X' = X — (, Z' = Z — £ momentum fractions. 
We show that writing the kernels W^(X, Z) in terms of the fractions X, X' , Z, Z' in the region £ < {X, Z} < 1 
gives the functions W(X, X'; Z, Z') which have the symmetry property with respect to the interchange of initial 
and final partons: W(X, X'; Z, Z') = W{X',X;Z',Z). For ( = one has X = X',Z = Z' and the kernels 
Wq = o(X, Z) acquire the DGLAP form. In Sec. VI, we discuss the QCD evolution of the nonforward distributions. 
Qualitatively, the evolution can be described in the following way. Due to the DGLAP-type evolution in the X > ( 
region, the momenta of partons decrease, and distributions shift into the regions of smaller X. However, when the 
parton momentum degrades to values smaller than the momentum transfer r — £p, the further evolution is like 
that for a distribution amplitude: it tends to make the distribution symmetric (or antisymmetric) with respect to 
the central point X = (/2 of the (0, £) segment. In section VII, we briefly discuss two basic uses of nonforward 
distributions: deeply virtual Compton scattering and hard elastic meson electroproduction. In particular, we show 
how to combine the definition of the gluon distribution through the matrix element of the gauge-invariant gluonic 
operator G^ Q ,(0)i? a fc(0, z; A)G^ v {z) with the usual Feynman rules formulated for the vector potential A^. In Sec. 
VIII, we discuss possible sources of PQCD factorization breaking for hard elastic electroproduction processes, due to 
singularities at the end-points of the integration region. In particular, we emphasize the importance of establishing 
the ^(0) = property for the nonforward distributions. In Sec. IX, we compare our notations, definitions 
and terminology with those used by other authors (off-forward parton distributions H(x,t;;t) introduced by X. Ji 
[ [To[ and nondiagonal distributions f(xi,X2) defined by Collins, Frankfurt and Strikman Jl^]). Sec. X contains 
concluding remarks. 

II. FORWARD AND NONFORWARD DISTRIBUTIONS IN SCALAR TOY MODEL 

A. Introductory remarks 

The parton distributions Fq{X; t; /i) parametrizing nonforward matrix elements (p'\O(0, z)\p) of composite two- 
body operators O(0, z) on the light cone z 2 — depend on four parameters. In addition to the "usual" parton 
variable X specifying the fraction Xp of the initial hadron momentum p carried by the active parton (more 
formally, X may be treated as the Fourier-conjugate parameter to {pz)), the functions !F^(X;t; ^i) also depend on 
the invariant momentum transfer t = (p' — p) 2 , the skewedness parameter £ = {rz)/(pz) (where r = p — p') and the 
evolution/factorization scale /i. The latter characterizes the subtraction procedure for singularities that appear on 
the light cone z 2 — (in general, fi may be different from the scale \ir introduced by the i?-operation for ordinary 
UV divergences, but the usual convention is to take \i = fin). Furthermore, depending on the type of the composite 
operator O(0, z), one would get quark, antiquark, flavor-singlet, flavor- nonsinglet, gluonic, spin-dependent, spin- 
independent, etc. distributions. In this situation, we propose to follow a step by step approach. We will start with 
simplest examples and then gradually proceed to more complicated ones. For this reason, we consider first a toy 
scalar model. The lowest nontrivial level corresponds to one loop Feynman diagrams. The relevant integrals are 
easily calculable, and their study provides useful information about the structure of the nonforward distributions, 
especially about their spectral properties, because the latter are insensitive to numerators of quark and gluon 
propagators and other complications brought in by the spin structure of the realistic QCD case. 

B. Forward distribution functions 

Our starting point is the scalar analogue of the usual "forward" parton distribution functions f(x). Consider a 
one-loop box diagram for a scalar version of the virtual forward Compton amplitude (Figs. la, b). Both incoming 
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and outgoing virtual "photons" have momentum q = q' — (p, where g' and p are lightlike momenta {q') 2 = 0,p 2 = 0. 
The "photons" couple with the constant e to a massive scalar "quark" field <f>. The initial and final hadrons are 
imitated by massless scalar particles with the momentum p. Their coupling to the quarks is specified by a constant 
g. In these notations, q 2 = —Q 2 = —2((pq'). Since (pq) = (pq') : the parameter £ coincides with the Bjorkcn 
variable £ — XBj = Q 2 /2(pq). Using the a-representation for the scalar propagators 



1 



£ ia(k -m +ie) da ^ 



m 2 — k 2 — it 

and calculating the resulting Gaussian integral over the loop momentum k we obtain for the diagram la: 



2„2 /■oo 



T a (p,q) = / <-xi>^- 



/ t\ «3 ~ C("2 + a 3 + a 4 ) 2 

2(pq)a 1 ■ ■ ■ A(m 

ai + a 2 + a 3 + «4 



da.lda2da.3dai 

A 2 



(2.2) 



We use the shorthand notation A = a\ + a 2 + 0:3 + a±. The large-Q 2 asymptotics is determined by integration 
over the region where the coefficient accompanying 2(pg') vanishes. Otherwise, the integrand rapidly oscillates 
and the result of integration is exponentially suppressed. Integration over a\ ~ region is evidently the simplest 
possibility. Other variants are a\ + a 2 + a 3 + a± ~ or 013 — C( a 2 + ot 3 + oti) ~ 0. It is easy to check that the 
leading power behavior is generated by the a\ ~ integration, which gives 



Ta{p, q) = -j/j / - - 7= ~ da 2 da 3 da 4 + 0(1/Q 4 ) , (2.3) 

16tt 2 J 2(pq')(a 3 /\-( + ie) X 2 



where A = 012 + ^3 + a 4- Introducing the distribution function 

i g 2 f°° f a 3 \ e -&(™*-ie) 

f( x ) = TT~2 S ( x ! — ~~ i da 2 da 3 da A , (2.4) 

^ Jo \ a 2 + a 3 + a A ) A 2 

we can write the leading power contribution in the parton form: 

TTiP, q) = - f ol ru a 1 ■ \ fit) dx = - f f , f{x) dx = f 1 t a (xp, q) f{x) dx . (2.5) 
Jo 2{pq ){x - C + »e) Jo ixp + q) 2 +te 







At the last step, we introduced the parton subprocess amplitude 



e 2 

t a (xp,q) = — . . (2.6) 

(xp + q) z + le 

Hence, the parameter x can be treated as the fraction of the initial momentum p carried by the quark interacting 
with the virtual photon. Note that the limits < x < 1 necessary for this interpretation of x are automatically 
imposed by the a-representation of fix). A similar result holds for the u-channel diagram 16: 

Tb*(p,q) = / r . ( C * . f(x)dx = - f ^—-fix)dx^[ t b (xp,q)f(x)dx. (2.7) 

Jo 2 \Pq)ix + Q-it) Jo ixp-qY+ie J 



The distribution function fix) is defined here by the same a-parameter integral (2.4). The latter can be easily 
calculated to give 

m = TA-2^-^0(0<x<l). (2.8) 



Note, that fix) is purely real. Due to singularity at x — £ in Eq. ( J2.5| ), the total amplitude T =T a + Tb has both 
real and imaginary parts. Since x > and ( > 0, its imaginary part is given by the s-channel contribution T a (p, g) 
only: 

h Imras ^' *> = I! Im ^ *> dx = if 2(^7) 6{x - f{x) dx =^§) = 2W)iS^^-^- (2 - 9) 
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The real part of T is given by TJ, and by the real part of T a : 

ReT™(p,q) = f Ret a (xp,q)f(x)dx = --^—p[ dx , (2.10) 

Jo 2 (PQ) Jo X ~C 

where P stands for the principal value prescription. 

To translate these results into the OPE language, we write the contribution of the diagrams la, 6 in the coordinate 
representation: 

Tip, q) = ( (p\ 0(O)0(«) \p) (e-' i(9z) + e^>) D m iz 2 ) d 4 z . (2.11) 



The large-Q 2 asymptotics of T(p, q) is given by the leading light-cone behavior of both the quark propagator 
Dmiz 2 ) = l/4i7r 2 (z 2 — ie) + . . . and the matrix element (p\ </>(O)0(z) \p) 

T( P ,q) = J ^(J-T^ WWW) \P)y-odh +Q(1/Q 4 ). (2.12) 
Defining the parton distribution function f(x) by 

(p | 0(O)0(«) | p) U 2 =o = J \ (e~ ix + e»^) fix) dx , (2.13) 



we rederive the parton formulas ( |2.5| ), ( |2.7| ). Basically, the integral ( 2.1 3| ) can be treated as a Fourier representation 



for the light-cone matrix element (p | <j>(Q)4>(z) | p) | z 2 =0 = fipz) which is a function of the only variable (pz). 
However, to derive the spectral constraint —1 < x < 1 for the Fourier partner of (pz) and establish the property 
fix) = f(—x), one should incorporate the fact that fipz) is given by Feynman integrals with specific analytic 
properties and that we have the same scalar field 4> at both points and z. The a-representation which we used 
above is one of the most effective (though perturbative) ways to take these properties into account. In Ref. pS}] 
(see also Sec. Ill below), the a-representation was used to prove that the constraint < x < 1 in Eq. ( p. 13 ) and 



similar (but more complicated) constraints for multiparton distributions and distribution amplitudes hold for any 
Feynman diagram. Two other approaches to studying spectral properties of parton distributions are described in 
refs. HQ. 

Anticipating comparison with the nonforward distributions discussed below, it is worth emphasizing here that 
the Bjorken ^-parameter is not present in Eq. ( |2.13j ) defining the parton distribution function fix). It appears 
only after one calculates the Compton amplitude Tip, q). 



C. Double distributions 



Now, consider a one-loop box diagram for the scalar analogue of the deeply virtual Compton scattering amplitude 
(Figs.lc, d). Using the same basic light-cone momenta p, q' as in the forward case, we write the momentum of the 
incoming virtual photon as q = q'—t^p. The outgoing real photon carries the lightlike momentum q' . The momentum 
conservation requires that the final hadron has the momentum (1 — £)p, i.e., in this kinematics we have a lightlike 
momentum transfer r = Qp. Since the initial momenta q,p are identical to those of the forward amplitude, the 
parameter £ coincides with the Bjorken variable XBj = Q 2 /2(po). In the a-representation, the contribution of the 
diagram lc is 

The large-Q 2 limit is again governed by the small-ai integration which gives 

ie 2 a 2 r°° 1 -i\(m 2 -ic) 
T c ip,q,q') = - 1 ^ - - - , -- - da 2 da 3 da 4 + 0(l/Q 4 ). (2.15) 



r,/ / \ «3 ~ C(«3 + (H) 2 

2(pq ) ai A (m — ie) 



a.\ + a-i + a 3 + a4 



da i da.2 da 3 da a 
A 2 



(2.14) 



G 



In the forward case, the ratio a^/X was substituted by the variable x which was interpreted then as the fraction of 
the initial hadron momentum carried by the active quark. The result expressed by Eq. ( 2.15| ) contains also another 
ratio uij\. So, let us introduce the double distribution 



F(x,y) 



16tt 2 



a 3 



Q'2 + 0.3 + Q '4 



a 2 



a 2 + a 3 + a 4 



-i\(m 2 — ie) 



A 2 



■ da2da^da4 . 



(2.16) 



It is easy to see that both variables x,y vary between and 1. Furthermore, their sum is also confined within 
these limits: < x + y < 1. Hence, F(x, y) = 6(x + y < 1) F(x, y). Using Eq. ( 2.16 ), we write the leading power 
contribution of T c (p, q, q') in terms of the double distribution: 



q,q') = - [ I 
Jo Jo 



—, — 7 Fix, y) dx dy 



■ F[x, y) dx dy = 



lo Jo ( X P + yr + q) 2 + ie 
The parton subprocess amplitude t c is given by 



l r i 



o Jo 



t c (xp + yr, q, q') F(x, y) 6(x + y < 1) da: dy . 



(2.17) 



t c (xp + yr,q,q') = - 



(xp + yr + q) 2 + ie 



(2.18) 



Hence, the momentum xp + yr of the quark interacting with the virtual photon originates both from the initial 
hadron momentum p (term xp) and the momentum transfer (term yr) . In a similar way, for the u-channel diagram 
Id, we get 



i r i 



i or /w — i — 7 — ^ F ( x > dx d y 

o Jo 2{pq'){x + yC, - ie) 



— F(x, y) dx dy 



lo Jo {xp + yr~ q') 2 + ie 
with the same double distribution F(x, y) given by Eq. ( 2.1 6| ) . In the explicit form, 



t d (xp + yr, q, q') F(x, y) 6(x + y < 1) dx dy , (2.19) 



F{x,y) 



16ir 2 m 2 



>{0<x + y<l) 



(2.20) 



Again, F(x,y) is purely real. Comparing the a-representations for f(x) and F(x,y), we obtain the reduction 
formula for the double distribution F(x,y): 



l-x 



F{x,y) dy = f(x) 



(2.21) 



Due to the restrictions x > 0,y > 0, the imaginary part of the total amplitude T = T c + Td is given by the 
s-channel contribution alone: 



-LlmT c ((,Q 2 ) ' 



2(pq) 



^0 



S(x + y( - ()F(x, y) 6(x + y < 1) dx dy 



1 



2C(P9) Jo 



F(x, 1 - x/Q dx = ^ J F(y(, y) dy = ^-r 1>(C). 



(2.22) 



The last form is similar to the expression for ImT in the forward case: one should just use the function 3>(C) instead 
of /(C)- Moreover, the integral defining $(£) looks similar to that appearing in the reduction formula ( 2.21 ). Still, 
the two integrals are not identical and, in general, $(C) ^ /(C)- Using the explicit form of F(x,y) for our toy 
model, we obtain 



*(C) 



9 

167r 2 m' 



■#(0 <C< !)• 



(2.23) 
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The factor (1 — Q present in f(Q (see Eq. (2.8) ), does not appear here. Note, however, that the difference is 
small for small £. 

In the OPE language, the basic change compared to the forward case is that we should deal now with the 
asymmetric matrix element (p — r \ <j)(Q)(j)(z) \p). Our definition of the double distribution F(x,y) corresponds to 
the following parametrization 

(p- r | <KO)0(z) |p)Ua=o = / / \(e- lx{pz) ^ v{rz) +e lx ^- i y {rz AF{x lV )9{x + y<l)dxdy. (2.24) 

Jo Jo 2 V / 

Taking the limit r = in Eq. ( |2.24| ) gives the matrix element defining the usual parton distribution function 



f(x), and we reobtain the reduction formula (2.21). Again, this definition of F(x,y) can be treated as a Fourier 
representation for a function of two independent variables (pz) and (rz), with the spectral constraints x > 0, y > 0, 
x + y <1 dictated by the analytic structure of the relevant Feynman integrals. An important feature implied by the 



representation (2.24) is the absence of the £-dependence in the double distribution F(x, y). The asymmetric matrix 
element 

(HI, of course, has C-dependence. But it appears only through the ratio (rz)/(pz) of variables in the 
exponential factor. In this treatment, £ characterizes the "skewedness" or "longitudinal momentum asymmetry" 
of the matrix elements. The fact that for the deeply virtual Compton amplitude T the parameter £ coincides 
with the Bjorken variable xsj = Q 2 /2(pq) is a specific feature of a particular process. The matrix element itself 
accumulates a process-independent information and, hence, has quite a general nature. 

Thus, despite the fact that the momenta p and r are proportional to each other r = Qp, there is a clear distinction 
between them, since p and r specify the momentum flow in two different channels. For r = 0, the momentum flows 
only in the s-channel and the total momentum entering into the composite operator vertex is zero. In this case, 
the matrix element coincides with the usual distribution function. The partons entering the composite vertex then 
carry the fractions x\ (i = 1,2) of the initial proton momentum. In general, — 1 < X4 < 1, but when X4 is negative, 
we should interpret the parton as going out of the composite vertex and returning to the final hadron. In other 
words, Xi can be redefined to secure that the integral always runs over the segment < x < 1 . In this parton 
picture, the spectators take the remaining momentum (1 — x)p. On the other hand, if the total momentum flowing 
through the composite vertex is r, the matrix element has the structure of the distribution amplitude in which the 
momentum r splits into the fractions yr and (1 — y)r = yr carried by the quark fields attached to that vertex. In 
a combined situation, when both p and r are nonzero, the initial quark has momentum xp + yr, while the final 
one carries the momentum xp — yr. Both the initial active quark and the spectator carry positive fractions of the 
lightlike momentum p: x + £y for the active quark and x — (y — (l—x — y) + (l — Qy for the spectator. However, 
the total fraction of the initial momentum p carried by the quark returning the fraction xp into the hadron matrix 
element is given by x — y( and it may take both positive and negative values. 



D. Asymmetric distribution functions 



Since (rz) — ((pz), the variable y appears in eq.( 2.24 ) only in the x + y( = X combination, where X can 



be treated as the total fraction of the initial hadron momentum p carried by the active quark. Since £ < 1 and 
x + y < 1, the variable X satisfies a natural constraint < X < 1. Integrating the double distribution F(X — y(,y) 
over y gives the asymmetric distribution function 

r x/( f x/c 
T C [X) = 9(X > () / F(X - y(, y) dy + 9(X < () / F(X - y(, y) dy , (2.25) 
Jo Jo 

where ( = 1 — £. The basic distinction between the double distribution F(x, y) and the asymmetric distribution 
function J-^(X) is that the former is a universal function in the sense that it does not depend on the skewedness 
parameter £ while the latter is explicitly labelled by it. Hence, we deal now with a family of asymmetric dis- 
tribution functions J r c i (X) whose shape changes when £ is changed. In our toy model, ^(X) has the following 
a- representation: 
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ig 2 f°° ( a 3 + C"2 \ e - lA (™ 2 - le ) 



T C {X) = SIX f » =- da 2 da 3 da A . (2.26) 

167r 2 7 V a 2 + a 3 + a4/ A 2 



Taking the integrals, we get the explicit form 



• FjY, = 16^{f () (0<X<O+ 1 T ^r 0kl A' 1^. ,2.27, 



One can see that when £ — > 0, the limiting curve for J-^(X) reproduces the usual distribution function: 

T C =o (X) = f(X) . (2.28) 



In general, this formula also follows directly from the definition of ^(X) and the reduction formula (2.21) for the 
double distribution F(x,y). 

The fraction (X — £) = X' of the original hadron momentum p carried by the "returning" parton differs from 
X by £: X — X' = £ [Q. Since X changes from to 1 and £ 7^ 0,1, the fraction X' can be either positive or 
negative, i.e., the asymmetric distribution function has two components corresponding to the regions 1 > X > Q and 
< X < £. In the region X > £ (Fig. 2a), where the initial parton momentum Xp is larger than the momentum 
transfer r = Qp, the function Tq{X) can be treated as a generalization of the usual distribution function f(x) for 
the asymmetric case when the final hadron momentum p' differs by Qp from the initial momentum p. In this case, 
Tq{X) describes a parton going out of the hadron with a positive fraction Xp of the original hadron momentum 
and then coming back into the hadron with a changed (but still positive) fraction (X — Q)p. The parameter £ 
specifies the longitudinal momentum asymmetry of the matrix element. 

In the region X < ( (Fig. 2b), the "returning" parton has a negative fraction (X — £) of the light-cone momentum 
p. Hence, it is more appropriate to treat it as a parton going out of the hadron and propagating along with the 
original parton. Writing X as X = YQ, we see that both partons carry now positive fractions YC,p = Yr and 
Yr = (1 — Y) r of the momentum transfer r. The asymmetric distribution function in the region X — YQ < £ looks 
like a distribution amplitude ^^(Y) for a 0^-state with the total momentum r = C,p: 

*cW= / F((Y-y)C,y)dy. (2.29) 
Jo 

In our model, 

yd Y ) = 7Tir^Y9(Q<Y <l). (2.30) 
Both F(x,y) and Tq(X) in our model are purely real. This result is determined, in fact, by general properties 



of the definition of the double distribution, Eq. (2.16). Indeed, let us introduce the Feynman parameters Pi by 
di = \f3i. After integrating over A, the only possible source of imaginary contributions is the denominator factor 
l/(m 2 — ie). However, since m 2 > 0, this factor is always positive and the integral is purely real. This would 
not happen if the initial "hadron" has a sufficiently large mass M > 2m. In this case, instead of l/(m 2 — ie) 
we would get l/(— M 2 (3^(l — /3 3 ) + rn 2 — ie) if the final hadron has the same mass M. Then, the denominator is 
not positive-definite if M > 2m, and the integral has both real and imaginary part. Clearly, the imaginary part 
appears because the initial hadron can decay into its constituents. If such a possibility is excluded, the double 
distributions F(x,y) and, hence, the asymmetric distribution functions !Fq(X) are purely real. 

The usual parton distributions f(x) are often related to imaginary parts, or more precisely, s- and w-channel 
discontinuities of parton-hadron amplitudes||. Note, that in our approach, the parton distributions are defined 
by form- factor-type matrix elements which depend only on momentum invariants p 2 , p' 2 , r 2 irrelevant to such 



'In a recent paper, L. Frankfurt et al. |23| discuss also discontinuities in the context of the nondiagonal distribution 
functions. 
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discontinuities (so far we even were treating these invariants as vanishing). The variable X in our definition only 
reflects a more complicated structure of the operator vertex. To illustrate this point, we write J-q(X) in the 
momentum representation (see Fig. 3a) 



5(X-(fcg')/W))<^ 



(27r) 4 J (fc 2 — to 2 + ie)((p — fc) 2 — m 2 + ie)((r — fc) 2 — to 2 + ie) 



(2.31) 



The S(X — (fcg')/(pg'))-function here corresponds to composite operator (denoted by a blob on Fig. 3a). Using the 
a-representation, one can take the Gaussian fc-integral and obtain the representation given by Eq. ( |Q6| ), which 



finally gives our purely real result ( |2.27|) 



It is worth emphasizing that the parton representations (2.5), (2.17) and (2.34) below are valid for the total 
Compton amplitude: there is no need to split the latter into its real and imaginary parts in order to define the 
parton distribution. To make a parallel with the traditional approach in which the parton distributions are defined 
through the discontinuities of parton-hadron amplitudes, let us calculate the fc-integral above using the Sudakov 
decomposition 



k = Zp + W + fe-L i 2(pq') = s 



(2.32) 



which gives 



2(2vr) 4 



d 2 fc. 



dij 



[Xrjs - k\ 



ie}[(X-l)r,s- 



k\ — to 2 



■ie]{(X-C)r,s-ki- 



(2.33) 



Looking at the location of singularities for the 77-integral, we immediately see that a nonzero result is obtained 
only when < X < 1. Furthermore, in the region ( < X < 1, the integral over rj is given by residue at 
7] = — (fc^ +m 2 —ie)/(l — X)s, which corresponds to substituting the ordinary propagator — l/[(p— fc) 2 — to 2 +ie] 
by the 5((p — fc) 2 — m 2 )-function for the line with momentum (p — fc). In other words, for £ < X < 1, our one- 
loop model for the function ^(X) is totally given by the residue corresponding to the s-channel cut through the 
parton-hadron scattering amplitude (see Fig. 36). On the other hand, in the region < X < the integral over 77 
is given by residue at 77 = (fc^ + m 2 — ie)/Xs, which corresponds to cutting the line with momentum fc (see Fig. 
3c). Such a cut cannot be related to s- or u-channel discontinuities In both cases, one can say that Tq(X) 
originates from a parton-hadron scattering amplitude T = iT whose imaginary part is given by one or another 
type of discontinuities. In our treatment, the only important fact is that the amplitude T is purely imaginary so 
that the distribution function Fq(X) is purely real. As we have seen above, the function J-^(X) can be written in 
several different ways, e.g., in the a-representation which can be integrated without taking any residues. 
In terms of T^{X), the virtual Compton amplitude T c +d{p, 9, q') can be written as 



T™ d (p,q,q') = 



2(pq) 



1 



X - C + ie X 



T Q (X)dX. 



(2.34) 



For a real function Tq(X), the imaginary part of Tjf? d (p, q, q') is determined by that of the short-distance amplitude 
(terms in square brackets). Since T^(X) linearly vanishes as X — > 0, the singularity 1/(X — ie) of the w-channel 
diagram Id gives a vanishing imaginary part. As a result, the imaginary part of the whole amplitude is generated 
by the 1/(A — £ + ie) singularity coming from the s-channel diagram lc: 



\ Im T C (C, Q 2 ) = tJ—, 
ne z 2{pq) 



8{X -QF^dX 



2(P9) 



mi 



(2.35) 



** I am grateful to L. Frankfurt for attracting my attention to this point and correspondence. 
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Hence, the integral <&(£) in Eq. (|2.22 ) is equal to .Ff (Q, i.e., to the asymmetric distribution function !Fq{X) taken 



at the point X = (. The parameter £ is present in ^(C) twice: first as the parameter specifying the longitudinal 
momentum asymmetry of the matrix element and then as the momentum fraction at which the imaginary part 
appears. As one may expect, it appears for X = XBj — C, j us t like in the forward case. Note, however, that the 
momentum (X — Qp of the "returning" parton vanishes when X = In other words, the imaginary part appears 
in a highly asymmetric configuration in which the fraction of the original hadron momentum carried by the second 
parton vanishes. Hence, ^(C) in general differs from the function /(C)- The latter corresponds to a symmetric 
configuration in which the final parton has momentum equal to that of the initial one. As discussed earlier, in our 
toy model /(C)/F"c(C) — f (0/^(0 = 1 — i- e -> -^cCO 1S larger than /(C), though the difference is small for small 
values of C- 

The fact that T^X) vanishes for X — has a rather general nature. Note, that for small X the function T^X) 
is given by its X < £ component 

r x/C 

T C (X)\ X < C = F(X-yC,y)dy. (2.36) 

Jo 

The size of the integration region is proportional to X and, as a result, T$ (X) vanishes like constx X or faster for 
any double distribution F(x,y) which is finite for small x and y. 

In the coordinate representation, the asymmetric distribution function can be defined through the matrix element 

(p' | </>(O)0(z) |p>U»=o = f \ (f iX ^ + e'^-OW) T^X) dX, (2.37) 
with £ = 1 — (p'z)/{pz). To re-obtain the relation between !F^{X) and the double distribution function F(x,y), 



one should combine this definition with Eq. ( 2.24 ). The £ — > reduction formula ( 2.28 ) trivially follows from Eq 



(2.37) 



Using translation invariance, we can write representation for a more general light-cone operator: 

(p' | 4>{uz)4>{vz)\p)\ z 2 =0 = J 1 (e-iXv(pz)+i(X-Ou(pz) + e -iXn(pz)+i(X-Ov(pz)'j jr^ X )dX. (2.38) 

This formula explicitly shows that if the parton corresponding to <p(vz) has momentum Xp, then the momentum 
of the parton related to <p(uz) is {X — C,)p and vice versa. 



E. Nonforward distributions 



Writing the momentum of the virtual photon as q = q' — is equivalent to using the Sudakov decomposition in 
the light-cone "plus" (p) and "minus" (q') components in a situation when there is no transverse momentum. An 
essential advantage of expressing the amplitudes in the a-representation is that it explicitly shows the dependence 
of the diagram on the relevant momentum invariants. This means that we can derive the parton picture both for 
zero and nonzero invariant momentum transfers t = (p' — p) 2 without bothering about an optimal choice of the 



basic vectors for the external momenta. Maintaining for simplicity p 



V 



0, we get for the diagram lc 



T c (p,q,q') 



eV 
16tt 2 



exp 



a 1 (2(pq)a 3 - Q 2 (a 3 + a 4 )) + ia 2 "4 . , 2 • \ 

A (m — it) 



Oil + «2 + «3 + a i 



daida2da^da4 
X 2 



(2.39) 



The small-ai integration then gives 



ie 2 g 2 



daida^da^ 



16^ 2 Jo 2{pq)(a 3 /X - ((1 - a 2 /X) + ie) A 2 



0{l/Q A ) , 



(2.40) 



where £ = Q 2 /2(pq) = xbj- Hence, introducing the ^-dependent double distribution 
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F(x, y; t) = ^ J™ S (x - a 3 /\) 5 (y - o 2 /a) 



^ita^a^jX— i\(m —ie) 



Hi"- Vn v ' I V ' 1 A 2 
we obtain the same parton formula, but with a modified parton distribution F(x, y; t) 



T c as (p,g,g') 



Jo Jo 



F(x, y; t) dx < 



(2.41) 



(2.42) 



2(pq){x + yC-( + ie) 

Moreover, the dependence on t appears only through the ^-dependence of F(x, y; t). Similarly, we can write down 
the a-representation for the u-channel diagram Id: 

ai(—2(pq')a 3 — Q 2 (a 2 + a 4 )) + to 2 «4 



eV r °° 



cxp < i 



A (m 2 — ie) 



da\doL2da-},da 4 
A 2 



Using 2(pq') = 2(pq) + t and integrating over small a\ gives the parton formula 



Tr( P ,q,q')= f f 
Jo Jo 



lo Jo 2(P<l)(x + y( + xt/2(pq)-ie) 
with the same i-dependent function F(x,y;t). In our model, 



F(x,y;t)dxdy 



F{x,y;t) = 



9(0<x + y<l) 



(2.43) 



(2.44) 



(2.45) 



16-7T 2 to 2 — ty(l — x — y) ' 

The parton subprocess amplitude in this case has the 0(t/2(pq)) = 0((t/Q 2 ) correction term which can be 
neglected in the large-Q 2 , fixed-i limit. Then the parton amplitude again depends only on the combination x + y(, 
and it makes sense to introduce the nonforward distribution 

r min{X/C,X/C} 

F(X-y(,y;t)dy, (2.46) 



T c (X;t) = f 
Jo 



which can be treated as the finite-i generalization of the asymmetric distribution function T^{X) (or more precisely, 
F${X) is the t = idealization of (X;t)). In our simple model, it can be calculated analytically: 



T c (X;t) 



g 2 At9(X>C) 



16tt 2 -txVTT^ 



In (t + Vl + r 2 ) + 9 (X < C) J ^ 



dy 



m 2 -ty(X- y() 



(2.47) 



where r = y/(-t/4m 2 )(l - X)/(l - (). The function T c {X;t) falls off with increasing \t\ like a form factor. 

The i-dependent distributions F(x,y;t) and !F^{X;t) in our model are purely real. Indeed, introducing again 
the Feynman parameters /3j by = A/3, and integrating over A gives the denominator factor l/(— 1/?2/?4 + ™ 2 — * e )- 
However, since t < 0, this factor is always positive and the integral is purely real. Imaginary part for F(x, y; t) and 
T^{X]t) would appear only if the initial hadron mass satisfies M 2 > 4m 2 . 

For real distributions, the imaginary part of the total Compton amplitude can be calculated by taking the 
imaginary part of the short-distance amplitude which picks out the function T^{C,\ t) 

„2 



(2.48) 



where T = y/(-t/4m 2 ){l - Q. 

In the OPE approach, the nonforward distribution is given by the matrix element 



(p'\Hmz)\p)i 



e tX(pz) + e i(X-0(pz) 



T c (X;t) dX. 



Taking the local limit z = 0, we obtain the following sum rule for T^{X;t) 



f 

Jo 



T c (X-t) dX = {p' | 0(0)0(0) | p) = F(t) , 



(2.49) 



(2.50) 



where F(t) is the toy model analogue of a hadronic form factor. 
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F. Time-like photon in the final state 



To give an example in which the skewedness parameter £ does not coincide with the Bjorken parameter XBj, 
let us discuss a situation considered in refs. [j]~T|— |l3|| , when the initial photon is spacelike q\ = q' — CiP while the 
final photon is timelike q 2 = q' + £2?- Here q' is a basic lightlike vector defining the Sudakov decomposition rather 
than the momentum of the final photon. Now, the Bjorken ratio given by XBj — — ?i/2(p5i) coincides with £1. 
However, by momentum conservation, the longitudinal part of the final hadron momentum differs from that of the 
initial one by (£1 + £2)^, i-e., the ^-parameter is given in this case by £ = Ci + £2, with £ 2 — q 2 /2{pq\). As a result, 
the scalar analogue of the timclikc photon elcctroproduction amplitude can be written as 

(for simplicity, we took here t = 0). Its imaginary part is proportional to the sum ^i+i^iCi) + •7 r c; 1 +£ 2 (£2)- Hence, 
having information about the imaginary part of such an amplitude for different values of £1 and £2, one can 
directly "measure" the asymmetric distribution function Tq{X) in the region X < £, where T$ (X) is similar to a 
distribution amplitude. 



1 



1 



X - £1 + it X-( 2 



T Cl+C2 (X)dX 



(2.51) 



III. ALL-ORDER ANALYSIS 

A. Handbag diagram to all orders 

Using the a-representation, one can write down the contribution of any diagram in terms of functions of the 
a-parameters specified by the structure of the diagram. Since the object of our interest is the matrix element of 
a two-body operator, we can extract it from the simplest handbag diagrams, i.e., those in which the q vertex is 
connected to the q' vertex by a single propagator. The contribution of any diagram of this type can be written as 
(see, e.g., MM) 



J. 2 a 1 A L {a) . a t A a (a) . a x A u (a) A t (a) ^ 2 1 
exp < iq + is + iu + it—— - 1 > ot a {m a -%()), (3.1) 

I LJ[a) LJ\ a ) t)(pi) LJ\ a ) I 

where s = (p + q) 2 , u = (p — q') 2 and t = (p — p') 2 are the Mandelstam variables, d is the space-time dimension, 
P(c.c.) is the relevant product of the coupling constants, z is the number of loops of the diagram and I is the 
number of its internal lines. Finally, D 7 A s , A u , A t , Al are functions of the a-parameters uniquely determined for 
each diagram. 

To describe them, we need definitions of a tree and a 2-tree of a graph. A tree (2-tree) of a graph G is a 
subgraph of G which consists of one (two) connected components each of which has no loops. Any tree G\ (2- 
tree G l 2 ) of G is determined by the set of lines a which should be removed from the initial graph G to produce 
G\ (G l 2 ). The product of the fto-parameters associated with these lines will be referred to as a-tree (a- 2- tree). 
The function D(a) is called the determinant of the graph. It is given by the sum of all a-trees of the graph G. 
By B{i\, . . . , i m \ji, ■ ■ ■ > jn) w e denote the sum of all a-2-trees possessing the property that the vertices i\, . . . , i m 
belong to one component, ,j n to the other, while the vertices not indicated explicitly may belong to either 

component. In these notations, 

a 1 A L (a) = B(q\p,q',p'); ai A s {a) = B{q,p\q' ,p'); a x A u {a) = B(q,p'\p,q'); A t {a) = B{q,q'\p,p') . (3.2) 

The mnemonics is straightforward: the square of the total momentum entering into one of the components (due to 
momentum conservation, it does not matter which one) just gives the relevant momentum invariant (see Fig. 4). To 
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get all the 2-trees corresponding to this invariant, one should make all possible cuts resulting in such a separation 
of external momenta. Note, that a.\ must be present in all terms of B(q\p, q',p'), B(q,p\q' ,p') and B(q,p'\p,q') 
because the vertices q, q' in these cases belong to different components. On the other hand, for B(q, q'\p,p') these 
vertices are in the same component. As a result, there are terms in A t (a) which do not contain a± as a factor, 
i.e., A t (a) = aiAf\a) + A[°\a) with A$(a) ^ and Af\a) ^ for a x = 0. Similarly, the function D(a) can 
be written as D(a) = ct\D\(a) + Do(a), where Di(a) is the determinant for the graph G\ obtained from G by 
deleting the line <7i, while D (a) is that for the graph Go resulting from G by contracting the line <7i into a point 
(and gluing the vertices q, q' into a single point). One can see that the function Do (a) can also be written in terms 
of the same a-2-trccs: 



D {a) = {B(q\p, q',p') + B(q,p\q',p') + B(q,p'\p, q') + B(q,p,p'\q')}/ai 

= A L (a) + A s (a) + A u (a) + A R (a) , 

where Ar(o.) is the function corresponding to the cut separating out the momentum invariant q a . 
leading large-Q 2 asymptotics, we integrate over the region ol\ ~ 0. This gives 



(3.3) 
To get the 



T as (p,q,q' 



diagr. 



P(c.c) 

(ATTi) zd / 2 



Y\ da a D d ' 2 (a) 



x exp < it 



4°V) 

D (a) 



a=2 



,A L {a) A s (a) A^ja) A?\a) 
D (a) ' D (a) D (a) D (a) 



ie) 



(3.4) 



(7 = 2 



Using s = -Q 2 + 2(pq),u = -2(pq) 
factor, we transform it into 



t and neglecting t and m 2 compared to 0(Q 2 ) terms in the denominator 



Q 



,A L {a) + A s {a) 



2(p<z) 



A„ 



A u {a) 



D (a) v ^ D (a) 

This expression has the structure similar to that of the one-loop contributions ( |2.40 ), ( 2.43| ). In particular, it can 



be converted into the form of the s-channel term ( 2.40 ) if we denote [-A s (a) — A u (a)]/Do(a) by x and [Al{o) 
A s (a)] /Do (a) by 1 — y. Analogously, to make it look like the u-channel term ( 2.43] ), we should take [-A s (a) — 
A u (at)\/ Dq(ol) — —x and [Ai,{a) + A s (a)]/Do(a) = y. If we want to have positive x, we should perform the first 
identification in the region where A s (a) > A u (a) and use the second one in the region where A s (a) < A u {a). In 
other words, we define the t-dependent double distribution by 



F{x,y;t) 



diagr. 



s i-i P{c.c.) 



(4iri) zd / 2 J 



Jo t-o D 0\ a ) 



- i^a a {m a - ie) 

o=2 



5 (1-y- 



A L (a)+A s (a) 



y - 



D (a) 
A L (a) + A s (a) 
D (a) 



A 3 (a)-A u (a) 
D (a) 
A u (a)-A s (a) 



D (a) 



6(A s (a) > A u (a)) 
9{A s {a) < A u (a)) 



(3.5) 



An intuitive interpretation is that when A s (a) > A u (a), the quark takes the momentum xp from the initial 
hadron. Its total momentum is xp + yr. Alternatively, when A s (a) < A u (a), the quark returns the momentum 
xp to the final state, and its total returning momentum is xp — (1 — y)r. Due to Eq. (3.3), we automatically 
have 0<x<1,0<j/<1. Furthermore, since x + y = [Al(o) + A u (a)]/Do(a) < 1 in the first region and 
x + y — [Ar(o) + A s (a)]/D (a) < 1 in the second one, we always have the restriction x + y < 1. 

Again, introducing the Feynman parameters Pi — oti/X and the common scale A given by the sum of all Ui- 
parameters, we can integrate over A to see that the resulting denominator factor 1/ (—tAf \a) / Do(a) + m 2 ) is 
positive for t < 0, and the double distribution is purely real. 

The same definition of the parameters x, y based on the a-representation can be used in the realistic case of 
spin-i quarks. However, one should take into account that the quark lines in that case are oriented. Depending 
on their direction, we should interpret the parton with momentum xp + yr either as a quark or as an antiquark. 
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The nonforward distributions J-^(X;t) can be obtained from the double distributions using Eq. flUf ). The 
restrictions x,y > 0, x + y < 1 guarantee that the total fraction X satisfies the basic parton constraint < X < 1. 
Furthermore, if the double distribution F(x,y;t) is finite for all relevant x,y, the nonforward distribution F^(X;t) 
vanishes (at least linearly) as X — > 0. 



(3.6) 



B. Alpha-representation and factorization 

Using the a-representation, we can write each perturbative diagram contributing to the virtual Compton scat- 
tering amplitude T(p, q, q') in any field theory model, including QCD (see Fig. 5) 

Pic c \ r°° -. 

TU(p,q,q') = i l ^-J^ J \[da a D- d l 2 {a)G{a,p,q,q';m <T ) 

/ a B L (a) + B s (a) B s (a) - B u (a) 
eXp \~ lQ D{a) + DW) 

f . B t (a)+B u (a) . 1 . r2 B 1 (a) + B 2 (a) v-^ . 2 . ,1 
exp |, f L_ + ;M ^ - , £ - *) ] 

The only difference is the presence of the preexponential factor G(a,p, q, q'; m a ) due to the numerator structure of 
the QCD propagators and vertices. It has a polynomial dependence on the momentum invariants. The functions 
B(a) are defined by the relevant 2-trees, e.g., B^(a) — B(q\p,q' ,p'), etc. 

In the region where Q 2 and 2(pq) — C.Q 2 are large, all the contributions having a powerlike behavior on Q 2 
can only come from the integration region inside which all the ratios Al/D, A b /D, A u /D vanish: if any of them is 
larger than some constant p, the integrand rapidly oscillates and the resulting contribution from such an integration 
region is exponentially suppressed. 

Since Al,A s ,A u and D are given by sums of products of nonnegative a-parameters, there are two basic pos- 
sibilities to arrange A^/D = 0. In the first case, called the "short-distance regime", A vanishes faster than D 
when some of the a-parameters tend to zero (small a correspond to large virtualities fc 2 , i.e., to "short" distances). 
The second possibility called the "infrared regime" , occurs if D goes to infinity faster than Aj when some of the 
a-parameters tend to infinity (large a correspond to small momenta fc, i.e., to the infrared limit). One can also 
imagine a combined regime, when Aj/D = because some a-parameters vanish and some are infinite. 

There exists a simple rule using which one can easily find the lines a whose a-parameters must be set to zero 
and those whose a-parameters must be taken infinite to assure that A^/D — 0. First, one should realize that 
Ai/D = means that the corresponding diagram of a scalar theory (in which G = 1) has no dependence on the 
relevant momentum invariant (Q 2 , s or u in our case). As the second step, one should incorporate the well-known 
analogy between the Feynman diagrams and electric circuits [ p8| : the a,j-parameters may be interpreted as the 
resistances of the corresponding lines a. In other words, a CT = corresponds to short-circuiting the line a while 
a CT = oo corresponds to its removal from the diagram. Hence, the problem is to find the sets of lines {<x}sD) {^Iir 
whose contraction into point (for {<t}sd) or removal from the diagram (for {ct}ir) produces the diagram which in 
a scalar theory does not depend on p 2 . Thus, the rule determining possible types of the powerlike contributions 
is the following: after the part of the diagram corresponding to a short-distance subprocess is contracted into a 
point and the part corresponding to soft exchange is removed from the diagram, the resulting diagram ("reduced 
diagram" , cf . |39],[l8| ) should have no dependence on large momentum invariants. 

Some examples are shown in Fig. 6. The simplest possibility is to contract into point some subraph H containing 
the photon vertices q, q' (Fig. 6a). The reduced diagram depends only on small invariants t, M 2 and masses m. The 
long-distance part corresponds to a nonforward distribution. This is the standard OPE configuration. However, 
since q' 2 is not a large momentum invariant: q' 2 — 0, there is a less trivial possibility shown in Fig. 66. In this case, 
there are two long-distance parts: one is given by a nonforward distribution again and the other can be interpreted 
as the distribution amplitude (hadronic component) of a real photon. Exchange of soft quanta between the two 
long-distance parts of Fig. 66 corresponds to a combined SD-IR regime (Fig. 6c): the a-parameters of lines inside 
H vanish while those belonging to the soft subgraph S tend to infinity. 
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One can easily invent other, more complicated configurations. Fortunately, not all of them arc equally important: 
different configurations have different Q 2 -behavior. The power counting is based on the observation that in the 
essential region of integration a a ~ 1/Q 2 for lines in the short-distance subgraph H and a a ~ Q 2 /p 4 for lines 
in the soft subgraph S (p 2 is some generic small scale, say, M 2 or to 2 ). In the momentum representation, this 
corresponds to k ~ Q for the -ff-lines and k ~ p 2 /Q for the Splines. As a result, in a theory with dimensionless 
coupling constants, we can use the dimensional analysis to derive that the contribution due to H behaves like 
Q-i-dn ^ ^ere du is the sum of dimensions (in mass units) of the fields associated with the external lines of H . We 
should also take into account extra numerator factors brought by these external lines. For instance, each external 
quark line adds a Dirac spinor u{p), two of them give u(p)u(p) ~ p, and p can combine with q from H to give 
(pq) ~ Q 2 . This means that each external quark line can bring an extra Q 1//2 factor. Note, that 1/2 is the spin of 
the quark. Similarly, an external gluon line can add a p^ factor. Combined with q^ from H it gives (pq) ~ Q 2 , 
i.e., the gluon line can bring an extra Q = Q 1 factor for the whole amplitude. Again, "1" is the spin of the gluon. 
Hence, each external quark or gluon line can give the factor Q Si ~ di = Q^ ti where U = dj, — s, is its twist. Note 
also that calculating the virtual Compton ampitude we do not convolute the vector indices jjL, v of the initial and 
final photon lines with momentum-dependent vectors. Hence, each external photon line gives only the factor Q~ l 
due to its dimension. Thus, the counting rule for the contribution of the hard subgraph H is 

t H (Q) < Q A - N -^ U , (3.7) 

where N is the number of external photon lines of the hard subgraph and summation is over quark and gluon 
external lines of H . For the simplest hard subgraph with two external quark lines this gives tjj(Q) ^ Q°, a scaling 
behavior as expected. For the configuration 66, the estimate is tjj(Q) ^ Q -1 - Hence, the contribution of Fig. 66 
is power-suppressed compared to that of Fig. 6a. Note that since the gluons have zero twist, the hard subgraph 
can have an arbitrary number of extra gluon lines without changing its power behavior. A similar power counting 
estimate ^2) based on k ~ p 2 /Q can be obtained for the soft subgraph S: 

ts(Q) < , (3.8) 

where the summation is over the external lines of S. Hence, exchanging a soft quark (Fig. 6c?) produces the l/Q 2 
suppression (S has then two external quark lines each having t — 1), while the exchange of any number of soft 
gluons is not necessarily accompanied by a suppression factor, at least on diagram by diagram level (for more 
details, see discussion in the next subsection). For the combined SD-IR configuration, the power counting estimate 
is 

t HS (Q) <Q 4 - n -^ hU Q-^ . (3.9) 
It is convenient to describe the power-low behavior of T(Q 2 ) in terms of the Mellin transformation 

ioo j 

T{c ^ = h! (10 *< J ) dJ - < 3 - 10 > 

— ioo 

Then the statement that T(Q 2 ) ~ (l/Q 2 )™ is equivalent to saying that the Mellin transform $(</) has a pole at 
J = —n. Take as an example the Mellin transform of the scalar diagram shown in Fig. la (it is essentially identical 
to the diagram lc): 

eV r 

(3.11) 

Small-ai integration corresponds to the simplest SD-regime 6a and generates the pole 1/(J+ 1) corresponding to 
the l/Q 2 asymptotic behavior. The relevant reduced graph is shown in Fig. 76. 



10L\- 



(a 3 + q 4 ) - a 3 /C 

Oil + 012 + 0(3 + 0:4 



exp {i t OL20i^l A — iX (to 2 — ie)} 



daida2da 3 da^ 



(011 



Ol2 + Qf3 



a 4 ) d/2 
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Another possibility to kill the dependence on large variables is to take «3 = a± = which corresponds to the 
reduced graph shown in Fig. 7c. To describe a simultaneous vanishing of two a-parameters, we use the common 
scale p = as + a± and the Feynman parameters 7^ — oti/p. The resulting p-integral p J pdp gives the pole 1/(J + 2) 
corresponding to a nonleading behavior 1/Q 4 . 

Furthermore, contracting the whole diagram into point (i.e. taking a, = for all a-parameters) we also obtain 
a reduced graph which does not depend on large variables. In this case, we introduce the common scale A = 
Oil + a 2 + c*3 + c*4 and the relative parameters Pi = a.i/\. In d = 4 dimensions, the integrand behaves like 
X J X 3 dX/X 2 which produces the pole 1/(J + 2) generating a nonleading behavior 1/Q 4 . However, if we take a 
scalar model in d = 6 space-time dimensions, then the integrand behaves like X J X 3 dX/X 3 and small-A integration 
generates the leading pole 1/(J+1). Note that in this case after the A-integration we still have the factor (3( capable 
of producing another 1/(J + 1) pole due to small-/3i integration. Hence, the total singularity of this diagram in 6 
dimensions is l/( J+ l) 2 , which gives T(Q 2 ) ~ (In Q 2 )/Q 2 . This corresponds, of course, to the scaling violation i.e., 
to evolution of the nonforward distribution. One can even extract the relevant evolution kernel from the remaining 
integral over p 2 , P3, Pi = 1 — Pi — Pz (the result, in fact, can be read off Eq. ( |2.27j ) ). Another observation is that 
if we simply integrate over small-cei region, the remaining integral da.1da.3da.il 'A 3 logarithmically diverges in the 
region of small A = a 2 + a^ + a^. This is the standard UV divergence of a matrix element of a light-cone operator 
in a theory with dimensionless coupling constants. 

Taking a 2 — > 00, we incorporate the IR regime corresponding to the reduced graph Id. If the quark corresponding 
to the (72 line is massless, the a 2 integral in this limit is a.^ 3 da 2 j a\. It produces the 1/(J+ 1) pole corresponding 
to the leading 1/Q 2 behavior. In the previous section, we did not see this contribution because the quark masses 
were assumed to be nonzero for all the lines. For nonzero mass, the factor exp[— ia 2 m 2 ] suppresses the large-a2 
integration and no poles in the J-plane are produced. In other words, the IR regime should be taken into account 
only for massless (or nearly massless) fields. Note, that in QCD the IR regime for the virtual Compton amplitude 
also gives 1/Q 2 behavior for massless quarks (see Eq. ( |3.8| ) ). However, in QCD this is a nonleading contribution 
compared to the scaling behavior produced by the purely SD regime 7b. 



C. QCD and gauge invariance 



After the SD- dominance is established, the next step is to write the contribution of the SD configuration in the 
coordinate representation (Fig. 8a) 

T(p, q, q') = f e-^ d 4 z [ (p' I <f>(z 2 )C(z, «i, js 3 M*i) I p) d 4 z 2 (3.12) 



(where is a generic notation for the quark fields -0, ip and the gluon field ^4) and expand the bilocal matrix 
element (p' \ 4>( z 2) ■ ■ ■ 4>{zi) \ p) in powers of (z-i — z\) 2 . Since we already know from the a-representation analysis 
that the virtualities inside the SD-subgraph are 0(Q 2 ), extra powers of [z 2 — z\) 2 for simply dimensional reasons 
result in extra powers of 1/Q 2 , and the leading large-Q 2 behavior will be given by the lowest term of this expansion 
corresponding to the lowest-twist composite operator. Parametrizing the nonforward matrix elements of the light- 



cone operators by formulas analogous to Eq. ( 2.38 ) gives the parton formulas similar to Eq. ( 2.34 ). Of course, 
this is just a general idea how to obtain the QCD parton picture for the SD-dominatcd amplitudes. Its practical 
implementation depends on specific properties of a particular process under consideration. 

The most important complication in QCD is due to the gauge nature of the gluonic field. In the Feynman gauge, 
the gluon vector potential has zero twist, and we should perform an infinite summation over the external gluonic 
lines both for the SD-subgraphs H and infrared subgraphs S. Consider the sum of gluon insertions into the quark 
propagator. It is well-known (see, e.g., ^-^3|) that after summation 

S c (ti - tj) + J S c (i - z)^gA^z)S c {z -n)d 4 z + ... = E(£, m A) S c (£ - ry) [l + 0(G)] (3.13) 
all the A-fields can be accumulated in the path-ordered exponential 
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E^,r];A)=Pexp \ig J* A„(z)dz»j (3.14) 

while 0(G) term depends on the gluonic fields only through the tensor G^ v and its covariant derivatives. Since 
G^y is asymmetric with respect to the interchange of the indices /j,, v, it should be treated as a twist-1 field. For 
the simplest SD configuration possessing a single long-distance part, combining the E- factors of all internal lines 
of the SD-subgraph, one gets gauge-invariant operators, e.g., q(zi)j u E(zi, z^; A)q{z2). 

If the lowest-order SD-configuration contains two long-distance parts (like in Fig. 66) , the gluonic corrections in- 
clude insertions into the external lines of the SD subgraph 86. The resulting path-ordered exponentials E n (oo, £, ; ^4) 
then go to infinity along the relevant light-cone directions, e.g., q' or p in case of hard electroproduction processes. 
However, for color-singlet channels there are at least two such exponentials and their long-distance tails cancel 
each other so that only the factors -£(£,77; A) related to SD-subgraph vertices £,77 remain. The basic effect of the 
exponential factor E(^,r];A) is that expanding operators C(^,ry) into the Taylor series, e.g., 

00 

q((ih»m,mA)q(v) = Y,-i AU1AU2 ■■■ A "' 1 mivD Vl D V2 ...£>„„ <z(0 ; A = r, - £ (3.15) 

n=0 

one gets local operators q^/ u D ul D^ 2 . . . D Vn q containing covariant derivatives D v = d v — igA v rather than ordinary 
ones. 

The cancellation of E n (£, 00; A) factors is very important for the success of the standard factorization program. 
Only after such a cancellation, the long-distance information is accumulated in universal matrix elements of gauge 
invariant light-cone operators. To illustrate the difference between color-singlet and color- nonsinglet channels, 
consider matrix element J(p,q') = (0\E q >(po, 0; A)ip(0)\p) of the quark field i(j(0) coming out of a state with 
momentum p and taken together with the accompanying gluonic field A which is then absorbed by a q' channel 
quark collecting the gluonic yl-fields into the E q >(oo, 0; A) factor (see Fig. 8c). Note that the latter can be written 
as 

(/>oo \ />oo />00 f-t 

\ A(t)dt)=l+ A{t)dt+ A(t)dt A{t 1 )dt 1 + ... 
Jo J Jo Jo Jo 



= 1 + / A(t)E q ,(t,0)dt (3.16) 
Jo 

where A(t) = igq'^A^ftq'). Substituting this result into the matrix element 

/>oo 

J(p,g') = W(0)|p>+ / (0\A(t)E ql (t,0)tP(0)\p)dt, (3.17) 
Jo 

shifting the arguments of all fields in the second term by tq' and performing the Taylor expansion 

E q , (0, -t) ^(-tq') = ]T K —j r {q'D) n ^) (3.18) 


we can take the integral over t to get 

00 -n+l 

j( P , q') = (o]vb) - £ -—^mWDY^p) , (3.i9) 

n=0 ^ ' 

where all the fields are taken at the origin. In fact, since 

(fi 2 ) , (3.20) 

the right-hand side of Eq. ( |3.19|) does not depend on (pq r ) (cf. ifjjj ). Note, that the new representation for J(p, q'), 
unlike the original one, is not explicitly gauge invariant. However, the ip-tevm can be represented as 
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Wb> = <°l*||y^> ( 3 - 21 ) 

and we can combine it with the first term from the sum to get a term containing a covariant derivative D — d — igA: 

JM) - ~ E j^r^\ A ^DTm . (3.22) 

Repeating this trick, i.e., representing the term outside the sum as 
and combining this term with the n — 1 term from the sum one obtains 

/ / 7~) \ 2 00 ■ 71 -4- 1 

It is clear now that we can write J(p, q') in a manifestly gauge-invariant form (cf. [ffBf ): 

J ^') = & ^Mfeww 55 <°i (w)°°^ b) ' (3 ' 25) 

In perturbation theory, matrix elements (0\(q' D) n ip\p) for finite n have ultraviolet divergences which can be reg- 
ulated in a standard way, e.g., by the dimensional regularization. After renormalization, we get one-loop terms 
like g 2 7„ log /i 2 . However, the anomalous dimension j n contains the usual (Y]™ 1/j) term [p4f which behaves like 
logn for large n. Hence, taking the formal limit n — > oo one encounters a logarithmic singularity, which requires 
an additional regularization on top of dimensional regularization (cf. p6| ) . The parameter characterizing the extra 
regularization can be taken proportional to [i, i.e., matrix element (0|.Ey(oo, 0; A)ip(0)\p) is the simplest example 
of a long-distance object with a double- logarithmic dependence on the UV cutoff p7j. Such objects ("collinear" or 
"jet" factors ^,|^|,|^|) play an important role in PQCD studies of Sudakov effects. However, within the standard 
factorization approach, presence of noncancelling double logarithms of Q 2 (reflected by double logarithms log 2 /i 2 in 
long-distance matrix elements) is treated as a failure of the factorization program, since the amplitudes in that case 
cannot be written through a convolution with parton distributions defined through matrix elements of light-cone 
operators [which have a single- logarithmic dependence on /i] . 

Another signature of Sudakov effects is the presence of the IR contributions (see Fig. 6c). Again, since all the 
hadrons participating in a hard exclusive scattering process are color singlets, summing over all soft gluon insertions 
one would get a path-ordered exponential over a closed contour, and by Stokes theorem 

(0|Pexp|i 5 ^A p (z)^||0) = l + (0|O(G)|0), (3.26) 

where 0{G) depends on the gluon field only through the field strength tensor which has nonzero twist generating 
a power suppression of the net IR regime contribution. 



IV. NONFORWARD DISTRIBUTIONS IN QCD 
A. Quark distributions. 

Let us discuss now the nonforward parton distributions in the realistic QCD case. For quarks, we should take 
into account that the field ip a {z) contains both the a-quark annihilation operator and the a-antiquark creation 
operator, i.e., the matrix element of the same light-cone operator tp a (0) . . . ip a { z ) determines distribution functions 
both for the quark and antiquark. Another complication is related to spin. There are two leading-twist operators 
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(0)7^(0, z; A)ip a (z) and V'a(0)7 M 7s-E(0, z; A)tp a (z), where, as discussed above, E(0,z;A) is the path-ordered 
exponential ( |3. 14 ) which makes the operators gauge- invariant. In the forward case, the first one gives the spin- 



averaged distribution functions f(x) while the second one is related to the spin-dependent structure functions g\{x). 
In this paper, we will concentrate on the ip a "fp,E(0, Z] A)tp a operators and gluonic operators with which it mixes 
under evolution. The relevant nonforward matrix element can be written asp| 

(p', s' | M0)zE(0, z; A)M?) I V, * >U==o (4-1) 

dX 

x dX, 

where M is the nucleon mass and s, s' specify the nucleon polarization. Throughout the paper, we use the "hat" 



= u(p', S ')zu(p, s) j (e- iX ^^(X; t) - e^ x -^pz)T«{X- 1)) 



(rather than "slash") convention z = z M 7 M . In Eq. (4.1), the quark and antiquark contributions are explicitly 
separated (cf. |5(]]). The exponential e~ lX< <pz) associated with the functions T^iX^f) and K,^(X\t) indicates that 
the field ip a (z) corresponds to the a-quark taking the momentum Xp from the nucleon. When the momentum Xp 
is taken from the nucleon by an a-antiquark, the corresponding annihilation operator is in ip a (0), and the functions 
,F a (X;i) and K, a (X;t) are accompanied by the exponential e l( - x ~^ pz ) corresponding to the momentum at the 
-0 a (z)-vertex. The antiquark terms come with the minus sign because the creation and annihilation operators for 
them appear in the reversed order. 

As emphasized by X. Ji juj, the parametrization of this nonforward matrix element must include both the 
nonflip term described by the functions ^(X'.t) and the spin-flip termP] characterized by the functions which we 
denote by K.^ (X;t). Taking the O(z) term of the Taylor expansion gives the sum rules (see [[To|) 

f 1 [^(X;t) - T«(X;t)] dX = F?(t) , (4.2) 
Jo 

[ [JC a c (X; t) - K%{X; 1)] dX = F 2 a (t) (4.3) 
Jo 

relating the nonforward distributions !F^(X;t), K,^(X;t) to the a-flavor components of the Dirac and Pauli form 
factors: 

Y,e a F?(t)=F x (t) , ^e a i^(t)=F 2 (t), (4.4) 

a a 

respectively (see also []52| and [Q). The spin-flip terms disappear only if r = 0. In the weaker r 2 = t = limit, 
they survive, e.g., F^O) = k" is the a-flavor contribution to the anomalous magnetic moment. In the formal t = 
limit, the nonforward distributions !F^(X;t), IC^(X;t) reduce to the asymmetric distribution functions !F^(X), 
K/HX). It is worth mentioning here that for a massive target (nucleons in our case) there is a kinematic restriction 

m 

- t > ( 2 M 2 /(. (4.5) 

Hence, for fixed the formal limit t — is not physically reachable. However, many results (evolution equations 
being the most important example) obtained in the formal t = 0, M = limit are still applicable. 

In the region X > £, the initial quark momentum Xp is larger than the momentum transfer r — Qp, and we 
can treat T^{X) as a generalization of the usual distribution function f a {x). When (-> 0, the limiting curve for 
F((X) reproduces f a (X): 



"Two other definitions of the nonforward parton distributions in terms of matrix elements of composite operators proposed 
by X.Ji Jl^ and Collins, Frankfurt and Strikman are discussed in Sec. IX. 

The possibility of a spin-flip in nonforward matrix elements was discussed earlier in refs. |t>l|J52j] . 
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^ =0 (X)=/„(X) ; J? =0 (X) = f s (X). 



(4.6) 



The spin-flip asymmetric distribution functions JC^(X) do not necessarily vanish in the £ — > limit. However, 
the relevant nucleon matrix element u(p')(zf — fz)u{p) is proportional to £ and the spin-flip term is invisible in 
the forward case. 

In the region X < (, one can define Y — X/( and treat the function Tq{X) as a distribution amplitude ^"(Y). 
In particular, the nonflip part in this region can be written as 



(u{p')zu{p) / \e- lY{rz) T?{CY)-e- l(1 - Y)[ - rz) Tl{CY)\ dY = (u(p 



■)zu(p) / , 
Jo 



4^(Y) dY . 



where the distribution amplitude ^(Y) is defined by 



^(Y) = ^(YC)-^(YC). 



(4.7) 



(4.8) 



The function ^"(Y) gives the probability amplitude that the initial nucleon with momentum p is composed of the 
final nucleon with momentum (1 — Q)p = p — r and a gg-pair in which the total pair momentum r is shared in 
fractions Y and 1 — Y = Y. 



B. Gluon Distribution 



For gluons, the nonforward distribution can be defined through the matrix element 



(p' | ZnZuG* (0)E ab (Q, z; A)G b av {z) \ p) | z2=0 



(4.9) 



= u(p')zu(p) (z ■ p) 



u (p') Zr 21 J Z u(p)( z ■ P) 



-iX(pz) + el (X-C)(pz) 



F 9 c (X\t)dX 



,-iX(pz) + e i(X-0(p*) 



JC 9 {X;t)dX. 



The exponentials e~ lX ^ pzS> and e ltyX ~^^ pz ^ are accompanied here by the same function T^iX-^t) reflecting the fact 
that gluon and "antigluon" is the same thing. Again, the contribution from the region < X < £ can be written 
as 



(p')zu(p) (z-r) / e~ iY{rz) ^ 9 C {Y; t) dY + "/C" term, 



(4.10) 



with the generalized Y <-> Y symmetric distribution amplitude ^f(Y;t) given by 



* g c(Y;t) = -(F 9 JY{;t)+f g c (Y{;t) 



(4.11) 



In the formal t = limit, the nonforward distributions J- 9 (X; t), JC 9 (X; t) convert into the asymmetric distribution 
functions !F 9 (X), K, 9 (X). Finally, in the £ = limit, T 9 (X) reduces to the usual gluon density 



FL (x) = xf g (x). 



(4.12) 



C. Flavor-singlet and valence quark distributions 



In our original definition (4.1) of the quark distributions, the exponentials exp[— iX{pz)\ and exp[i(X — £)(pz)] 
are accompanied by different functions !Ff{X;t) and !F^(X;t)), respectively. In many cases, it is convenient to 
introduce the flavor-singlet quark operator 
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Oq(uz, vz) = (uz, vz) 



(4.13) 



where 



0^(uz,vz) = 



ip a (uz)zE(uz, vz] A)?p a (vz) — ip a (v z) zE(v z , uz; A)ipa{uz) 



(4.14) 



The nonforward distribution function ZF9(X; t) for the flavor-singlet quark combination ( |1.13| ) 



(p',s'\0 Q (uz,vz)\p, s)\ z 2 =0 = u(p',s')zu(p,s) 



i 

o 2 



^ — ivX (pz)-\-iuX' (pz) (pz)-iuX(pz) 



(where X' = X — £) can be expressed as the sum of "a + a" distributions: 

T^(X; t) = ]T(^(X; t) + T«(X; t)) . 



ff(X;t)dX + "/C" 
(4.15) 



(4.16) 



Writing the contribution from the < X < £ region as 

(u(p')zu(p)(z-r) [ e~ tY{rz) ^(Y;t)dY + term, (4.17) 
Jo 

we introduce the flavor-singlet quark distribution amplitude ^^(F;t) which has the antisymmetry property 
$®{Y; t) = — 9®(Y; t) with respect to the Y <-> Y transformation. 
Another combination of quark operators 



Ol\uz,vz) = \ 



ip a (uz)zE(uz, vz; A)ip a (vz) + ip a (vz)zE(vz, uz; A)ip a [uz) 



(4.18) 



corresponds to the valence combinations J-J a (X; t) = T"(X; t) — T"(X; t): 



{p' 1 s'\O { ~ ) {uz,Vz)\p,s)\ z 2= =u{p', S ')zu{p,s) I I e -ivX(pz)+iuX'(pz) + jvX>(p Z )- l uX{pz) Fj*(X,t)dX +«K. V 



C 

(4.19) 

In both cases (see Eqs.( l.lq ),( 4.19 ) ), two possible exponential factors are accompanied by the same distribution 



function, just like for the gluon distribution. In the region < X < £, the function T^ a (X;t) can be written 
in terms of the flavor-nonsinglet distribution amplitide ^Y a (Y;t) which is symmetric ^Y a (Y;t) = ^Y a (Y;t) with 



respect to the Y «-» Y interchange. 



V. EVOLUTION EQUATIONS FOR NONFORWARD DISTRIBUTIONS 

A. General formalism 

Near the light cone z 2 ~ 0, the bilocal operators (f)(0)<p(z) develop logarithmic singularities \nz 2 , so that the 
formal limit z 2 — ► is singular. Calculationally, these singularities manifest themselves as ultraviolet divergences 
for the light-cone operators. The divergences are removed by a subtraction prescription characterized by some 
scale /i: T^(X;t) — > T^(X;t; a). In QCD, the gluonic operator 

O g (uz,vz) = z^G^ a (uz)E ab (uz,vz;A)G b a „(vz) (5.1) 
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mixes under renormalization with the flavor-singlet quark operator. At one loop (i.e., in the leading logarithm 
approximation) , the easiest way to get the evolution equations for nonforward distributions is to use the evolution 
equation |3^,|34| for the light-ray operators^. For the flavor-singlet case, it reads 



a H> Jo Jo l 



u, v)Ob{uz, vz) Q(u + v < 1) du dv , (5-2) 

where v = 1 — v and a,b — g, Q. For valence distributions, there is no mixing, and their evolution is generated by 
the QQ- kernel alone. Inserting Eq. (5.2) between chosen hadronic states and parametrizing the matrix elements 
by appropriate distributions, one can get the well-known evolution kernels such as DGLAP and BL-type kernels 
and also to calculate the nonforward kernels R ab (x, y; £, 77) and W® b (X, Z). The kernels R ab (x, y; £, 77) govern the 
evolution of the double distributions: 

1 1 /• 1 

H—F a (x,y;t;ti)= \ / £ R ab (x, y; £, 77) F b (£, 77; t; /1) 6(£ + r, < 1) d£ d V , (5.3) 

b 



d/i 



Jo 



where a and b denote g or Q. Another set of kernels W^ b (X, Z) dictates the evolution of the nonforward distributions 
and asymmetric distribution functions: 

d f 1 
"A 1 Jo b 

The evolution of the double distributions will be briefly discussed later in Sec. VI. Here we will discuss the 
structure of the W£ b (X, Z) kernels. Since the form of the equation is not affected by the ^-dependence, "f" will 
not be explicitly indicated in what follows. 

Before starting the actual calculations, one should take into account that the gluon distribution Ff{X) is ac- 
companied by the sum of two exponentials while the flavor singlet quark distribution !F9(X) with which it mixes 
is accompanied by the difference. This sign change is, in fact, compensated by the extra (pz) factor in the right- 
hand side of the gluon distribution definition. The set of evolution equations for T®(X), !Ff {X) can be obtained 
by substituting the definitions of the gluon ( |4.9| ) and quark (4.15) distributions into Eq. (|5.2| ) and performing 
the Fourier transformation with respect to the (pz)-variable. For this procedure, the (pz)-factor is equivalent to 
differentiation d/dX while \/(pz) results in an integration over X. Note, that both operations change the relative 
sign of the exponentials. Hence, it is convenient to introduce first the auxiliary kernels Al£ b (X, Z) which would 
appear in the absence of the (pz) mismatch. They are directly related by 

M£ b {X,Z)= [ ( B ab (u,v)6{X -uZ + v{Z - ())9{u + v < l)dudv (5.5) 
Jo Jo 

to the light-ray evolution kernels |53]]34|] , which we write here in the form given in Ref. jn| : 
B QQ {u,v) = ^C F (l + 6(u)[v/v)+ + 5(v)[u/u]+ - b(u)S(v] 

7T \ Z 



B gQ {u,v) = ^C F [2 + 5(u)5{v) 
BQ g (u, v) = —Nf (1 + 4uv — u — v) , 

7T 

B gg (u,v) = ^N c L(l + 3uv-u-v) + ^5(u)5(v) + L(u) ^ - 8(v) ' 



(5.6) 



+ {u <-> v} 



^This procedure was also used in a recent paper [|35| . I was informed by J. Blumlein that its authors agree with my results 
for the W((X, Y) kernels given below. 

***We prefer to use the kernels B a t(u,v) which have the symmetry property B a t(u,v) = B a i,(v,u) and are related to the 
K ab (u,v) kernels of Ref. @ by B ab (u,v) = -K ab (u,v). 
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The W-kernels are related to the M-kernels by 

rQQlv 7\ _ n rQQ 



W 99 (X, Z) = Mf {X, Z) , W^{X, Z) = M^ W {X, Z), (5.7) 
^W 9Q (X,Z) = -M 9Q (X,Z) , W^ 9 (X,Z) = -^M^ 9 (X,Z). (5.8) 



Hence, to get W 9 (X, Z) we should integrate M 9 (X, Z) with respect to X. The integration constant can be fixed 
from the requirement that W 9 ®{X,Z) vanishes for X > 1. Then 

W 9Q (X, Z) = j M 9Q (X, Z) dX . (5.9) 



Integrating the delta- function in eq.( |5.5| ) produces four different types of the 9- functions, each of which corre- 
sponds to a specific evolution regime for the nonforward distributions. In two extreme cases, when £ = or £ = 1, 
the evolution equation reduces to known DGLAP and BL-type equations, respectively. 

B. BL-type evolution kernels 

When £ = 1, the initial momentum coincides with the momentum transfer and Tq{X) reduces to a distribution 
amplitude whose evolution is governed by the BL-type kernels: 

W^i{X, Z) = V ab (X, Z). (5.10) 



Taking £ = 1 in Eq. (5.5) we obtain 

M^ 1 (X,Z) = U ab (X,Z)= [ [ B ab (u,v)8{X-uZ-v(l-Z))6(u + v<l)dudv. (5.11) 

Jo Jo 

Eliminating the S- function, one would observe that in the regions X < Z and X > Z the U ab (X,Z) kernels 



are given by different analytic expressions. However, from the representation (5.11) and the symmetry property 
B a b(u,v) — B ab (v,u) it follows that U ab (X,Z) = U ab (X, Z). Hence, it is sufficient to know the /7-kernels in the 
X < Z region only. The basic function Ug b (X, Z) = 9(X < Z) U ab (X, Z) can be calculated from 

Ug b (X,Z) = I f B ab (v-(X -v)/Z,v)dv. (5.12) 
z Jo 

The total kernel U ab (X, Z) then can be written as 

U ab {X, Z) = 9(X < Z) U£ b (X, Z) + 9{Z < X)U$ b (X, Z) . 
One can easily derive a table of B — > Uq conversion formulas for all the structures present in the i?-kernels: 



S(u)S(v)^S(Z-X) , l->£ , S(u)-^0 , 5(u)(~) 

Z v \v J 



0. 



u ( X\ 1 „, , u A ( X\ 2 1 



Z Z - X ' w u \Z Z-X' 



x ( X \ X 2 (\ X x\ . . 

Using Eqs.(|5.6|) and this table, we can get the BL-type kernels V ab (X, Z). Before doing this, we note that the BL- 
type kernels appear as a part of the nonforward kernel W" (X, Z) even in the general £ ^ 1,0 case. As explained 
earlier, if X is in the region X < £, then the function J-q(X) can be treated as a distribution amplitude (Y) with 



24 



Y = X/C- For this reason, when both X and Z are smaller than £, we would expect that the kernels (X, Z) 
must simply reduce to the BL-type evolution kernels V ab (X/C, Z/Q. Indeed, the relation (^5|) can be written as 



Mf{X,Z) = - f [ B ab (u,v 
s Jo Jo 



) 6 (X/C - uZ/C, - v(l - Z/Q) 9(u + v < 1) du dv . 



(5.14) 



Comparing this expression with the representation for the UQ b (X, Z) kernels, we conclude that, in the region where 

If 



X/C < 1 and Z/Q < 1, the kernels Mf(X, Z) are given by 



Mf(X,Z)\ < {X!Z} < c = -U ab {X/C, Z/Q 



(5.15) 



From the expressions connecting the W- and M-kernels, we obtain the following relations between the non- 

■ward evolution kernels W^ b (X, Z) in the regioi 
W£ b (X,Z)\ < {x ,z}<() and the BL-type kernels V ab (X,Z): 



forward evolution kernels W" b (X, Z) in the region < {X, Z} < C (let us denote them by Lf(X, Z) 



L Q Q Q (X, Z) = - V™ (X/C, Z/Q ; Lf(X, Z) = V gQ (X/C, Z/Q ; 
Lf°(X, Z) = ± V** (X/C, Z/Q ; Lf(X, Z) = ± V" (X/Q Z/Q . 



Explicit calculations based on Eqs.(5_ L 5)-(5 : 9), ( 5.10 ), ( 5.16 ) give 

1 



Y w (x, z) = —c f 



z \ z-x 



V Q9 (X, Z) 



as 



X 
~Z 



4(1 -X) 



>(X <Z) 
1 - 2X~ 



f {X -> X,Z -> Z}| 
(X < Z) - {X -> X, Z -> 2} 



V9 Q {X , Z ) = ^C F ^2-^e(X<Z) + 9(X>Z)} 



V"(X, Z) = — N c { 2 ^- ( 3 - 2X 



+6(X - Z) 



1-X 

z 

' Po 

2N r 



1 



Z-X 
1 dz 



I) 



1 



] (X < Z) + {X X,Z ^ Z}. 



Note, that the V gQ (X, Z) kernel can be represented as the sum 

v gQ (x,z) = ^-c F + ^c F Ui-^\e(x <z) - {x -> x,z -> z}| 



(5.16) 

(5.17) 

(5.18) 
(5.19) 

(5.20) 
(5.21) 



of a constant term and a kernel which is explicitly antisymmetric with respect to the {X — > X, Z — > Z} transfor- 
mation. In fact, the constant term does not contribute to evolution since the flavor-singlet distribution amplitude 
tyQ(Z) with which it is convoluted is antisymmetric tyQ(Z) = —^f^(Z). For the same reason, the convolution of 
V g Q(X, Z) with ^Q(Z) determining the evolution correction to !F^(X) behaves like X 2 for small X. 

Furthermore, the BL-type kernels also govern the evolution in the region corresponding to transitions from a 
fraction Z which is larger than C to a fraction X which is smaller than C- Indeed, using the ^-function to calculate 
the integral over u, we get 



i r x/c ( \ 

M?(X,Z)\ X < C < Z = - J B ab ([l-X/Z-v(l-C/Z)},v)dv, 



(5.22) 



which has the same analytic form ( 5.12 ) as the expression for Mf(X, Z) in the region X < Z < C- For QQ, gg and 
Qg kernels, this automatically means that W? (X, Z)\x<q<z is given by the same analytic expression as Lf b (X, Z) 
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for X < Z. Because of integration, to get Wj? Q (X, Z) one should also know (X, Z) in the region ( < X < Z. 
However, our explicit calculation confirms that W^(X,Z) in the transition region X < ( < Z is given by the 
same expression as L 3 ®(X, Z) for X < Z. 

Note, that the evolution jump through the critical fraction £ is irreversible: the 5- function in Eq. ( 5.14 ) requires 
that X/C, = v + (1 — u — v)Z/( or X < £ if Z < (. To put it in words, when the parton momentum degrades 
in the evolution process to values smaller than the momentum transfer £p = r, further evolution is like that for a 
distribution amplitude: the momentum can decrease or increase up to the r-value but cannot exceed this value. 



C. Region Z > <, X > ( 

Recall, that when X > £, the initial quark momentum Xp is larger than the momentum transfer r = (p, and we 
can treat the asymmetric distribution function !Fg(X) as a generalization of the usual distribution function f a (X) 
for a somewhat skewed kinematics. Hence, we can expect that evolution in the region £<X<1,£<Z<1 
is similar to that generated by the DGLAP equation. In particular, it has the basic property that the evolved 
fraction X is always smaller than the original fraction Z. The relevant kernels are given by 



M? b (X, Z)\ 



C<x<z<i 



B ab [ [1 - X/Z - v(l - C/Z)} 



dv . 



(5.23) 



Changing the integration variable to w = v(l — £/Z)/(l — X/Z) = v/(l — X'/Z'), we obtain the expression in 
which the arguments of the _B-kernels are treated in a more symmetric way 



Mf{X,Z)\ c 



<X<Z<1 



Z -X 
ZZ' 



B ab (w (1 - X/Z) ,w(l- X'/Z')) dw , 



(5.24) 



where X' = X — £ and Z' = Z — ( are the "returning" partners of the original fractions X, Z. Moreover, since 
Z-X = Z'-X', the kernels Mf{X, Z) are given by functions symmetric with respect to the interchange of X, Z 
with X' , Z' . This observation can be used to check the results of calculations. However, since we are dealing with 
the asymmetric situation X > X',Z > Z' , other practical applications of this symmetry are not evident at the 
moment. Again, we can easily obtain a table for transitions from the i3 a f,-kernels to the M ah -kernels for the region 
C < X < Z < 1: 



uv 



S(u)S(v) 

Z -X 
~* 6ZZ' 



S(Z-X) ; 1 



X 
~Z 



'z 1 

S(u 



z-x 

ZZ' 



v 
v 



-S(v) 



(«• 



Z-X 
2ZZ' 



5(v) 



1 



Z-X 



z-x 

2 



~Z' 
X' 



z + z> 



z 7 



(5.25) 



Introducing the notation P£ b (X, Z) 
the P-kernels-^ : 



W£ b (X, Z)\(_<x<z<i and using the formulas given above, we calculate 



Expressions for the nonforward generalization of the DGLAP evolution kernels (in different notations) were given in 
the review jl2| by L. Gribov, E. Levin and M. Ryskin. They discuss the generalized DGLAP kernels in the context of the 
electroproduction amplitude with a timelike photon (or Z°) in the final state. However, as the longitudinal momentum 
asymmetry parameter ( for their kernels they took the ratio ("2 = (§2./2{j>q\) involving only the invariant mass g| of the final 
photon. As we have seen in Sec. II F, the correct value for ( in this case is C = Ci + C2, where (1 is the usual Bjorken 
parameter £1 = —qf/2(pqi) 



2G 



P? Q {X, Z) 



O c 



r 1 1 + z 2 1 1 
- S(X z) J ^— dzj - -p qq (x/z) . 



P? 9 (X,Z) = ^N f ' 

C V ' ' 7T / ZZ> 



Pf{X, Z) 



—P Qg (X/Z). 



P° 9 (X, Z) 




(5.26) 
(5.27) 
(5.28) 



^P gg {X/Z). 



(5.29) 



As expected, the P£ b (X, Z) kernels have a symmetric form. 
P£ b (X, Z) are related to the DGLAP kernels in the C 
one should take into account that the asymmetric gluon distribution function Tf (X) reduces in the limit ( = to 



The formally divergent integrals over u and v provide here the usual "plus" -type regularization of the 1/(Z — X) 
singularities. The prescription following from Eqs. ( 5.24 ), ( |5.2E| ) is that combining the l/(Z — X) and 8(Z — X) 
terms into \T^{Z) — !F^(X)]/ (Z — X) in the convolution of Pq{X, Z) with J-^iZ) one should change u — > 1 — X/Z 
and v -> 1 - X'/Z'. 

The arrows indicate how the nonforward kernels 
limit when Z = Z' and X — X' . Deriving these relations, 

Xf g (X) rather than to f g (X). 

In the region Z > £, the evolution is one-sided: the evolved fraction X is smaller than Z. Furthermore, if Z < £ 
then also X < Z, i.e., distributions in the X > £ regions are not affected by the distributions in the X < £ regions. 
Hence, just like in the DGLAP case, information about the initial distribution in the Z > £ region is sufficient for 
calculating its evolution in this region. This situation may be contrasted with the evolution of distributions in the 
Z < £ regions: in that case one should know the asymmetric distribution functions in the whole domain < Z < 1. 

Qualitatively, the evolution in the X, Z > ( region proceeds just like in the DGLAP evolution: the distributions 
shift to smaller and smaller values of A. In the DGLAP case, the distributions approach the 6(x) form condensing 
at a single point x — 0. In the asymmetric case, the whole region Z < ( works like a "black hole" for the partons: 
after they end up there, they will never come back to the X > £ region. Inside the Z < £ region, the evolution 
is governed by the BL-equation transforming the ^^(Y) distribution amplitudes into their asymptotic forms like 
YY,YY(Y — Y) for the quarks and (YY) 2 , (YY) 2 (Y — Y) for the gluons; a particular form is dictated by the 
symmetry properties of the relevant operators. 



VI. ASYMPTOTIC SOLUTIONS OF EVOLUTION EQUATIONS 

A. Evolution of asymmetric distribution function 

To describe the qualitative features of the QCD evolution of the nonforward distributions, we will consider 
the simplest case, i.e., the evolution equation for the flavor- nonsinglet (valence) functions. Then Qg, gQ and gg 
kernels do not contribute, and the evolution is completely determined by the QQ-kernel. The multiplicatively 
renormalizable operators in this case were originally found in Ref. |^] 

o n = {zd + ) n i>\ a zc z j 2 (z d /zd+)i>. (6.i) 

Here we use the symbolic notation (z d / zd + ) of Ref. ||, with d = d — d,d+=d + d and Cn 2 {y) being the 
Gegenbauer polynomials. This means that the Gegenbauer moments 

C c (n, fi) = [ C 3 / 2 (2Z/( - 1)T C (Z; M ) dZ (6.2) 
Jo 

of the asymmetric distribution function T^{X\ /j) have a simple evolution: 



27 



C f (n,/i) = C ( (n,no) 



In /xq/A 
In /i/A 



(6.3) 



where /3 = 11— |iVj is the lowest coefficient of the QCD (5- function and 7„ are the nonsinglet anomalous dimensions 

HP 



1 



n+1 



(n + l)(n + 2) 



2 £j 



i=2 



(6.4) 



For n = 0, the Gegenbauer moment coincides with the ordinary one and, since 70 = 0, the area under the curve 
remains constant. Other Gegenbauer moments decrease as fj, increases. For the ordinary moments of the nonforward 
distribution 



M N (C,fx) 



f T C (X;^)X N dX , 
Jo 



(6.5) 



using explicit expression for the Gegenbauer polynomials we can derive the following expansion over the multi- 
plicatively renormalizable combinations C^(n, /i): 



N 



M N ((^)=C N AM(iY+l)!]T(-l)' 



2(2n + 3) 



n=0 



(N + n + 3)\{N - n)\ 



C f (n,/x) . 



(6.6) 



We can also write the expression which gives the evolved moments M.n{Q, A 4 ) m terms of the original ones: 



N 



M N ((, f i) = C N n\ (N + iy.J2 



(-l) n 2(2n + 3) 
^ (N + n + 3)UN- n y. 

n— 



In no/ A 



In n/A 



7n/ 'Pa n 



E 

k=0 



{-l) k (k + n + 2)\ 
2( k kl (k + 1)1 (n - k)l 



(6.7) 



With increasing N, the number of contributing Gegenbauer moments C^(n,fi) in Eq. fl6.6| ) increases. An 
important observation is that the nonevolving (and (^-independent, but i-dependent) term C(0) contributes to each 
moment. As a result, in the fi — > 00 limit, all the moments tend to constant values determined by the n = term 
in the sum fl6.6j): 



M N {C,,V -> 00) = C 



A' 



{N + 2)(N + 3) 



C(0) 



C(0) 

c 



6(A7C)(1 - x/ox" dX 



(6.8) 



Note, that the last integral involves only the X-values smaller than £. This means that in the limit /x — > 00, the 
function J-^(X; fi — > 00) completely disappears from the region X > i.e., it reduces to the distribution amplitude 
^(F) which ultimately tends to the usual asymptotic shape 6y(l — Y) in the Y = X/C, variable: 



T C (X; M - 00) = 6C(0)X(1 - X/O/C 2 



(6.9) 



One may also be interested in finding expressions showing how the function !Fq (X; fj,) changes its shape from an 
arbitrary original curve ^(X;^) to the asymptotic one. Note, that the Gegenbauer moments for £ < 1 involve 
integration regions in which the argument C^!" 1 '(2Z / '£ — 1) of the polynomials extends beyond the segment (—1, 1) 
where they form an orthogonal set of functions. Hence, a formal inversion of the Gegenbauer moments is only 
possible for ( = 1. In this case, the inversion produces the standard solution of the evolution equation for a 
distribution amplitude |^J| 



^ c=1 (Jf;/x) = ]r 



71 = 



4(2n + 3) 
(n + l)(ra + 2) 



xxc z J 2 {2x- i) 



In no/ A 




ln/i/A 


Jo 



Cl' 2 {2Z -l)T c=1 {Z-^)dZ . (6.10) 
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Thus, if the initial distribution coincides with one of the eigenfunctions XXCn (2X — 1), the evolution is very 
simple: the function just decreases in magnitude without changing its form. An attractive feature of such a 
situation is that approximating the initial distribution amplitude by a few lowest Gegenbauer polynomials one 
obtains a simple model of its evolution. Inspired by this observation, one may be tempted to construct a similar 
representation for the evolution of the asymmetric distribution function. Using the expansion of the light-cone 
operator ■ip(0)\ a z'ip(z) over the multiplicatively renormalizable operators O n (see |Q ) 

^(0)X a z^(z) = V"(_i)" 2 ( 2n + 3 ) [ (uu) n+1 O n (uz)du (6.11) 

^— ' 77.1 /„ 



n=0 "' 1/0 

and inserting it into the nonforward matrix element, we obtain 

^ c (JC;m) = £(-1)" 2 -^^{j^Y' P0 CCdn^o) j\uur^S^(X-uOdu. (6.12) 

Integrating (uu) n+1 S (n \X - u() over u, we get the Gegenbauer polynomials C^" 1 '{2X / Q — 1) accompanied by the 
spectral condition X < (. This means that the formal integration does not give a correct result for functions which 
do not vanish outside the region X < Q. For such functions, one should first perform the summation over n (which 
is, of course, practically impossible) and only then take the u-integral. 

Another limit in which the integral over u can be taken safely is £ = 0. For small £, the Gegenbauer polynomials 
are dominated by the senior power Z n and in the £ — > limit we obtain 

T^{X^) = Y j t^L 5 (n) {x) (^JAy n ' Pa J^T c=Q {Z-^)Z n dZ, (6.13) 

i.e., the well-known result that the moments of the usual parton densities have a simple DGLAP evolution. Note, 
that in this case, the functions which evolve without changing their shape are S^ n \x). From a pragmatic point 
of view, this observation is of little use. Modeling the solutions of the DGLAP equations is known to be a rather 
complicated excercise usually involving a numerical integration of the evolution equations. 



Hence, the representation (6.12) should be understood only in the sense of (mathematical) distributions in X 



rather than functions. To get meaningful results, one should integrate them over X with some smooth function. 



In particular, integrating it with X , one obtains the formula (6.6) for the evolution of the X moments of 
nonforward distributions. 



B. Evolution of double distribution 



Solving the evolution equation for the valence double distribution F (x, y; /i) defined by 



(p-r, S '|0<-)(O,z)Up, s)\ z2=0 



u(p',s')zu(p,s) 



e -ix(pz)-iy(rz) , e ix(pz)-iy(rz) 



F(x, y; n) 8(x + y < 1) dx dy , 



(6.14) 



we can give an alternative derivation of the asymptotic form of the valence nonforward distribution !F^[X; /i). The 
/i-dependence of F(x, y; fi) is governed by the evolution equation 



d f f 

M -T- F (x, y; m) = / / Rqq(x,v;C,v) f (^, 

"A* Jo Jo 



r]]n)dr). 



(6.15) 



Since the integration over y converts F(x,y) into the parton distribution function f(x), whose evolution is governed 
by the DGLAP equation 
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the kernel Rqq(x,u; £,,r]) must have the property 

" 1 

For a similar reason, integrating -Rqq(x, y; £, n) over x one should get the BL-type kernel: 



RQQ{x,y;i,r])dx = V QQ {y,r]). 
Explicit calculation gives for Rqq(x, y; £, 77) the following result 

V, C, »?) = ^C F i 1 0(0 < as/e < min^/n, y/f?}) - - z/0<% - «) 
9(0 < x/i < l)x/£ 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



-6(x/€ - y/rj) + ^=6(x/£ - y/rj) 



26(1 - x/i)8(y - rj) 



(WO 

It can also be obtained from the kernel Bqq(u, v) using the relation 

RQQ(x,y,^r)) = jB QQ (y - r)x/£,y - tjx/£). 



o 1 - a 



■ dz 



(6.20) 



It is easy to verify that the spectral constraint x + y < 1 is not violated by the evolution: the kernel Rqq(x, y; £, 77) 
has the property that :c + y<lif£ + r)<l. Using our expression for i?QQ(x, y; £, 77) and explicit forms of 
the Pqq(x/0 and Vqq (y, 77) kernels (see Eqs. ( 5.26 ), ( 5.17 ) ) one can check that RQQ(x,y; £,77) satisfies the 
reduction formulas ( |6.17| ) and fl6.18| ). To solve the evolution equation, we combine the standard methods used to 
find solutions of the underlying DGLAP and BL evolution equations. To solve the DGLAP equation, one should 
consider the moments with respect to x. Multiplying Eq. (6.15) by x n , integrating over x and utilizing the property 

Rqq(x,v;€,v) = RQQ( x /£,,y; hv)/(, w e get 



(If I 

where F n (y; fi) is the nth cc-moment of F(x, y; fi) 



d f 1 
H-7-F n (y;n) = / R n (y,r])F n (ri; fj,)dr) , 
an Jo 



Fn(y;n) 



1-3/ 



and the kernel R n (y,r)) is given by 



R n (y,v) = — C F 

7T 



{ 


/ \ n+1 


1 1 




.CO 


n + 1 77 — y 



x n F(x, y; fi)dx 



% <v) + {y^y,v^v} 



8(y-rj) 



dz 

z 



(6.21) 



(6.22) 



(6.23) 



It is straightforward to check that R n (y,r]) has the property 

Rn(y,v) w n(v) = Rn{n,y)w n {y), 

where w n (y) — (yy) n+1 ■ Hence, the eigenfunctions of R n (y,r]) are orthogonal with the weight w n (y) — (yy) n+1 , 
i.e., they are proportional to the Gegenbauer polynomials C^ +3 ^ 2 (y — y) (cf. Now, we can write the general 

solution of the evolution equation 



F n (y;n) = 2(yy) r 



„ +1 (2n + l)!(2n + 2)! ^ (2fc + 2n + 3) fc! cn+3/2 



E 



n\ (n+1)! ^ (2n + k + 2)\ 



(y - y) 



log(MoM) 
log( M /A) 



A nk (n ), (6.24) 
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where 



F n (y^ Q )C- +3/2 (y-y)dy 



(6.25) 



and the anomalous dimensions 7! are related to the eigenvalues of the kernel i? n (y, 77). They coincide with the 
standard nonsinglet anomalous dimensions 7jv (|6.4[ ): 7^ = Jn+k- Since 7q°' ) = 0, while all other anomalous 
dimensions are positive, in the formal /i — > 00 limit we have Fo(y,/j, — > 00) ~ yy and F n (y,n — > 00) = for all 
71 > 1. This means that 



F(x,y;fi^ 00) - 5(x) yy, 



(6.26) 



i.e., in each of its variables, the limiting function F(x,y;fx — ► 00) acquires the characteristic asymptotic form 
dictated by the nature of the variable: S(x) is specific for the distribution functions 55 5(|, while the yy-form is 
the asymptotic shape for the lowest-twist two-body distribution amplitudes S0] . For the asymmetric distribution 
function this gives ^(X, [i — > 00) ~ (X/( 2 )(l — X/Q. This result was already obtained in the previous subsection. 



VII. BASIC USES OF NONFORWARD DISTRIBUTIONS 



A. Deeply virtual Compton scatterring 



Using the parametrization for the matrix elements of the quark operator, we can easily write a parton-type 
representation for the handbag contribution to the DVCS amplitude: 



2(pq 



E e " [(-^ + j^wf +pV*)) { fi(p')5'«(p)Tg(o + ^(p'x^ - rq'Hpmo 

(7.1) 



where <j' = 7 M g' M , and T a (£) are the invariant amplitudes depending on the scaling variable £. In particular, 



X - C + ie X - ie 



(^(X;t)+^(X;t)) dX. 



(7.2) 



Since the nucleon is the lowest bound state in the 3-quark system, the nonforward distribution function for 
t < is real. Hence, the imaginary part of Tp(C) can be produced only by singularities of the terms in the square 
brackets. Taking into account that the nonforward distributions vanish for X = 0, we conclude that only the term 
containing l/(X — £ + ie) generates the imaginary part: 



1 



(7.3) 



with a similar expressions for ImT£ G P (Q. As discussed in Sec. I, the function .F^ (£; t) does not coincide with the 
usual parton distribution f a ((), even in the formal t — > limit. To get the real part of the 1/(X — £ + ie) terms, 
one should use the principal value prescription 



ReT£(£) = -P f\^(X;t) + ^(X;t)) 
Jo 



dX 

x~( 



(7.4) 



Since the principal value prescription is based on cancellation of X < £ and X > £ parts of the integral, it makes 
sense to preserve J^{X;t) as a single function. Splitting it into X < £ and X > £ components, one would simply 
get two divergent expressions for the real part of the amplitude. 
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Let us study how these formulas are modified by the evolution. At one loop, the lnQ 2 term can be easily 
calculated using the coordinate representation: 



Ti(p,q,q') 



2tt 



-i(qz) 



{p'\ip(uz) 



2iw(z 2 ) 2 



ip(vz) | p) lnz 2 Bqq(u, v) du dv . 



(7.5) 



Parametrizing the matrix element by the nonforward distribution (|4. 1| ) , we obtain for the s-channcl short-distance 
amplitude 



t[ s) (p,q,q') = In Q 2 



r , . Bqq(u,v) dudv 

yqp) 



l 



o (q + {vX - uX')p) 2 +ie 2 



= 7><2 



1 pi 



Bqq(u, v) du dv 



(1 - u)X' -vX + ie' 



where X' = X — £. Using explicit expression for the Bqq(u,v) kernel, we obtain 



1 



X' + ie 



*!(*) = "v7 



1 



X' + ie 



3 , 
2 +ln 



X' + ie 

-( + ie 



lnQ 5 



A similar expression can be derived for the evolution of the w-channel-type term: 



X-ie 



t\{X) = 



1 



X-ie 



1 + f^ 



2+ In WO 



lnQ 5 



(7.6) 



(7.7) 



(7.8) 



Clearly, the it-channel term can be obtained from the s-channel one by the change X' — > X, C — > — £. In the 
region A < £, both i" and if are real. Furthermore, it is easy to establish that the correction terms in both cases 
vanish when integrated with the asymptotic distribution 6A(1 — A/£)/£, explicitly showing that the latter does 
not evolve with Q 2 . Note that tf(X) is purely real in the whole range < X < 1, while £f(A) is purely real only 
in the region A < C- For X > £, the latter has an imaginary part: 



t\(X) = -P 



1 



x-c 



+ nr6(X -Q + -1Q 



2tt 



F S 77 



1 



x-c 



+ITT 



\o(x>CY 




rin|A/C-l|" 




[ (x-C) _ 


+ 


[ x-c 


.) 



+ mS(X - C) 



lnQ 2 



(7.9) 



This information can be used to write down the expression showing the leading logarithm evolution of the 
function (C', Q 2 ) determining the imaginary part of the amplitude: 



^c(C; Q 2 ) = ^c(C; Ql) + ^c F in q 2 /q 2 f { S(x - Q ( § - £ ^ 



x-c 



F c {X-Q 2 )dX. (7.10) 



Evidently, the expression in the braces is given by the nonforward evolution kernel (C, X) (5.2£). For the usual 
distribution function the analogous equation contains the DGLAP kernel P(C/X): 



f(C;Q 2 ) = f(C;Q 2 ) + ^C F inQ 2 /Q 2 J j 



6(x - 



- 2 



1 dz \ 1 + (C/A)- 
o l-z) X-C 



f(X;Q 2 )dX. (7.11) 



The comparison of the two expressions shows that evolution of the function ^(C',Q 2 ) is n °t identical to that of 
f(C; Q 2 ). Recall also that in the forward case the lowest-order amplitude is proportional to l/(X — C + *e) + 1/(A + 
C-ie). 



B. Gluonic contribution to hard exclusive meson electroproduction 

The kinematics of hard exclusive meson electroproduction processes -f*P — ¥ Mp' is very close to that of the virtual 
Compton scattering, especially in a situation when one can neglect the mass of the final meson. Again, one can use 
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the a- representation rules to determine possible regimes capable of producing a powerlike contribution for large Q 2 . 
The basic difference is the absence of the regime analogous to short-circuiting a subgraph containing the photon 
vertices, since instead of the final photon described by an elementary field we have now a bound state. Hence, 
the leading short-distance regime corresponds to contraction into point of a subgraph which contains the virtual 
photon vertex and located in the middle between the two long-distance-sensitive pp'- and g'-components of the 
diagram. The pp'-component is described by the nonforward distribution function while the g'-part is parametrized 
by the meson distibution amplitude. 

Depending on the type of lines connecting the short-distance subgraph with the (p' \ . . . \p) matrix element, one 
deals either with quark (Fig. 9a) or gluonic (Fig. 10) contributions to the lowest-order amplitude. The structure 
of the quark contribution is similar to that of the hard-gluon-exchange contribution to a meson electromagnetic 
form factor, with the distribution amplitude of the initial state substituted by the quark nonforward distribution. 
There is also an analogue of the soft contribution to the meson form factor (see Fig. 96). ft corresponds to the 
infrared regime a aa — > oo. 

Let us concentrate here on the gluonic contribution which requires a proper handling of restrictions imposed by 
gauge invariance. Using the coordinate representation for the hard propagators, we can write the contribution of 
Fig. 10a as 

T«(p,p',q') = J (q',M\m 7 a S c {-z 1 )T a ^S c {z 1 -z 2 )r b ^(z 2 )\0) (p'|A«(z 1 )^(z 2 )|p) d 4 Zl d 4 z 2 , (7.12) 

where r a ,r fc are the SU(3) color matrices. The first matrix element here can be expressed through the meson 
distribution amplitude <^(t) while the second one is related to the asymmetric gluon distribution. Other 3 lowest- 



order diagrams can be written in a similar way. Applying formally the power counting (see Eq. (3.7) and discussion 
preceding it) , we may conclude that each gluonic contribution has an extra Q 2 factor compared to the quark term, 
since the quarks have twist 1 while the twist of the gluon vector potential is zero. Technically, the enhancement 
appears when the p^Pv factor from the matrix element (p'\A'] 1 (zi)A h v (z2)\p) combines with the g^, q' u factors from 
hard propagators and polarization vectors, thus producing the estimate (p'|v4A|p) ~ Q 2 . However, the power 



counting formulas like (3/7) only give an upper estimate for the relevant contribution. The actual behavior is 
determined by the twist to of the composite operator O constructed from the elementary fields corresponding 
to the external lines of the SD-subgraph. It is well-known that the simplest gauge-invariant composite operator 
containing two gluonic fields is G iip G p vl and its twist equals 2 rather than 0, just like for the lowest-twist ?/> . . . ip 
operator. Diagrammatically, this means that, in the Feynman gauge, the leading power terms of four lowest-order 
diagrams completely cancel each other and the total result is suppressed by 1/Q 2 compared to leading contributions 
of separate diagrams. In general, picking out nonleading power terms (higher twist contributions) is a notoriously 
difficult problem of perturbative QCD. However, in our case, the cancellation of leading terms is guaranteed by 
gauge invariance of the total result. Hence, choosing a gauge in which the combination q'^q' v {p'\A a ^(zi)A b v (z 2 )\p) 
is prevented from producing the (q'p) 2 factor, we would eliminate the artificially enhanced terms on diagram by 
diagram basis. This goal is achieved if one uses the gauge q'^A^z; q') — 0. Then A^ can be expressed in terms of 
the field-strength tensor G MP (see, e.g., [5^]) 

A li {z;q') = q'P ( G^z + aq') e^ a da. (7.13) 



This representation also makes it easy to parametrize the matrix element (p'\A'^ l (zi)A b J (z2)\p) in terms of the 
gauge-invariant gluon distribution: 

/ /i ah, is Ah, /xi \ i S ab u(p')q'u(p) ( Pu,q' v +Pi,q' u \ .„ 

(p'\A^ Zl ;q')Al(z 2 ;q') \p)\( Zl -^ = 2(q> ■ p) {~ 9 ^ + (p ■ q') J {TU) 



X 



1 1 r 

_ e -iX(pzx)+iX (pz 2 ) , e iX (pz 1 )-iX(pz 2 ) 

o 2 . 



f a c {X) 

■dX + "AC". 



(X - ie)(X' + ie) 



In Ref . H , the amplitude of hard diffractive electroproduction was calculated for the longitudinal polarization 
of both the virtual photon (e^« = (g' M + Cp fJ ')/Q) an d produced vector meson (cy = g' M /my). In this case, the 
contribution of Fig. 10a in the (g'^4) = gauge can be written as 
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'•' ( CP - t£ (X-Qp + T? 

((p-T q i)^((X~()p + Tqiy 



T£ L (p, q', r) ~ u(p')q'u(p) / dr Vv (r) / Sp 



-fl^H 7 ^\ J 7y r^Ty — rr dA, (7.15) 

IP'S) / (X-ie)(X-C + ie) 

where w(t) is the distribution amplitude of the longitudinal vector meson. This gives 

u(p>)q'u(p) f 1 dr f 1 gffl 

<2™v Jo T Jo - C + *e) 

Other diagrams give similar contributions, differing only in the r-dependent factor. For Fig. 106, one should 
substitute 1/t by 1/f, while Figs. 10c, d both have 1/rf factor. Since 1/r + l/f = 1/rf, the total contribution also 
has the 1/rf structure 

T L L{p,q,r)~— / lp v {t)— " — — dX , (7.17) 

Qmy J tt J X{X-C + ie) 

where i/l — £ comes from u(p') = VI — C &(p)- The amplitude Tll(p, q,r) has imaginary part due to the factor 
l/(X-{ + ie): 

-ImT Li (C)~^^(C) r^rfr. (7.18) 



7T (QlTlV 







In Ref. the gluonic matrix element was approximated by the gluon distribution function f g (Q). To get our 
result from that of Ref. ||, one should substitute there f g (C) by \fl ~ (J 7 ^ (()/(. 

Though the asymmetric distribution function !F^{X) coincides with Xf g [X) in the limit £ = 0, in general these 
two functions differ when C 7^ 0. As discussed earlier, the imaginary part appears for X = £, i.e., in an asymmetric 
configuration in which the second gluon carries a vanishing fraction of the original hadron momentum, while Qf g (£) 
corresponds to a symmetric configuration in which the final gluon has the momentum equal to that of the initial 
one. 



VIII. FACTORIZATION AND END-POINT EFFECTS 
A. General remarks 

The standard question about PQCD applications for hard processes is whether factorization of short- and long- 
distance contributions is maintained in higher orders. Since the Feynman integrals can be written in different 
representations, one can approach the factorization problem in various ways. In particular, the classic studies 
of deep inelastic scattering in QCD |5^,|55||36| relied on the operator product expansion in which the coordinate 
representation plays a crucial role. The claims that factorization also holds for a more complicated Drell-Yan process 



60 1 were supported by studies |6l| , |39| , |62| based on the analysis in the momentum representation (see, however, 
[47 1). The early studies of exclusive processes in QCD which started with the analysis of the large-Q 2 behavior 
of the pion EM form factor also incorporated both the OPE-like coordinate representation methods [^|j4| and 
momentum- representation oriented approaches ||,|6| . Factorization was intensively studied in the following years 
(see [ fl9|]64] ] and references therein) . Referring an interested reader to Ref. |M for a recent momentum-representation 
analysis of factorization for hard exclusive electroproduction processes, here we briefly discuss possible sources of 
factorization breaking analysing them within our approach || based on the combined use of the a-representation 
and the OPE-type methods. 
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B. Structure of the lowest-order term 



Exclusive processes are rather vulnerable to factorization breaking. In contrast to inclusive cross sections, 
factorization for exclusive ampitudes may fail even at the tree level. Hence it is a good idea just to write down 



the lowest-order contribution and carefully look at it. Take the DVCS amplitude (7.2). It has terms 1/(X — ie) 
and 1/(X — £ + ie) which are singular for X = and X = respectively. An immediate question is whether 
these singularities appear within the region of integration and if yes, whether they are inside that region or at 
its end-points. To be prepared to address this question, we performed a detailed study of spectral properties of 
nonforward distributions. Our a-representation analysis shows that < X < 1. Since the singularity 1/(X — (,+ie) 
is inside the integration region, we can write it as P{1/(A — £)} — iir5(X — £): it generates both real and imaginary 
parts of the amplitude. On the other hand, the 1/{X — ie) singularity is at the end-point, and the relevant real 
part is given by a divergent integral unless the nonforward distribution J-^(X) vanishes at X = 0. Hence, to claim 
factorization for the real part, it is absolutely necessary to give the arguments that ^(0) = 0. In our analysis, 
we derived !F^{X) from the double distribution F(x,y). The basic expression for Fc,(X) shows that !F^{X) <~ X 
for any F{x,y) which is finite as x,y — > 0. One can get J^(0) ^ only if F{x,y) is singular for x = 0, e.g., if it 
behaves like 5(x) and does not vanish when y = 0. If F(x, y) has such a behavior, there should be a very special 
reason for it. 

Similarly, for the meson electroproduction, the integral over t contains the factor l/r(l — r) singular at the 
endpoints r = , r = 1. Again, the factorization formula makes sense only if the distribution amplitude <p(r) 
vanishes for r = 0, 1. Since <p(r) is analogous to the £ = 1 limit of a nonforward distribution, we may expect that 
it also vanishes at t = because of small phase space for the r — > configuration. Furthermore, since for massless 
quarks (p(l — r) = ±y>(r), if ip(r) vanishes at t — 0, it also vanishes for r = 1. 

Of course, even if the vanishing at end-points holds for any diagram of perturbation theory, this still does not 
mean that the nonperturbative functions have the same property. So, a cautious statement might be that if in 
perturbation theory some function does not vanish at a particular end-point, it is unlikely that it will vanish nonper- 
turbatively. If it vanishes perturbatively, there is some hope that this property is preserved for the nonperturbative 
function. 

The standard procedure to get an educated guess concerning the end-point behavior of hadron distribution 
amplitude tp(r, n) is to study the asymptotic fj, — > oo limit of their evolution. This idea is equivalent to saying that 
ip(r, fi) has the same behavior at the end-points as the relevant BL evolution kernel V(r, r'). In particular, in refs. 
J63[ it was shown that (^ as (r) ~ r(l — r) both for the longitudinally and transversely polarized p-mesons. 

Similar estimates of the end-point behavior of the distribution amplitudes follow from QCD sum rule consider- 
ations. In particular, if perturbative term W ert (r : M 2 ) of the QCD sum rule (M 2 is the SVZ-Borel parameter) 

f p (p(T)e~ m ^ M2 + higher states = W ert {r, M 2 ) + condensates (8.1) 

vanishes for r = and r = 1, one can argue that because of quark-hadron duality, tp(r) should also vanish at the 
end-points. For the correlators corresponding to the leading-twist p-meson distribution amplitudes, we have indeed 
W ert (r, M 2 ) - t, (1 - r) at the end-points. 

Note, furthermore, that both quark and gluon propagators of the simplest hard subraph have denominators 
proportional to r. However, for a longitudinally polarized virtual photon, only the 0(t) term in the numerator of 
the quark propagator survives which converts the 1/t 2 singularity of the hard amplitude into 1/t. This will not 
happen if the p-meson is transversely polarized. Hence, for transverse polarization one would face the integral with 
<Pt(t)/t 2 which logarithmically diverges if <pt(t) ~ t for small r. This result excludes the transverse case from 
straightforward PQCD applications. This fact was repeatedly emphasized in refs. P, |l^ , |l8| , [2"l| . 

C. Double-flow regime 

One of the lessons from the discussion above is that taking into account only the denominators of the "hard" 
quark and gluon propagators one is guaranteed to get a 1/t 2 factor capable of destroying factorization from the 
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very start. It is the cancellation of one power of r by a numerator factor in case of a longitudinally polarized virtual 
photon which makes the factorization possible. In the absence of this cancellation, e.g., for transversely polarized 
p-meson, even if we take <pt{t) ~ t(1 — r), the integral would logarithmically diverge. One may object that in 
such a situation factorization still works if <Pt(t) vanishes faster than r as r — * 0. Note, however, that evolution 
generates terms proportional to r: 



since V(t,t') ~ r for small r. In the presence of nonzero masses m or other infrared cutoffs, one should change 
1/r by 1/(t + m 2 /Q 2 ). As a result, the logarithmic divergence converts the r integral into an extra lnQ 2 /m 2 . 
Together with the evolution logarithm hiQ 2 /p 2 , they would amount to a double logarithm in a one loop diagram 
of Fig. 11 type. It should be emphasized that this is not a Sudakov double logarithm. In particular, in two loops 
one would only get In 3 Q 2 (In 2 Q 2 from evolution and InQ 2 from the r-integral) rather than In 4 Q 2 . The possibility 
to get an extra logarithm in the form-factor-type amplitudes was discovered a long time ago in a scalar model 
(see e.g., Ref. |6E{| ), In a scalar model, there are no numerator factors to moderate the 1/t 2 singularity, hence 
such a possibility is always realized. In Ref. ||, the diagram of Fig. 11 type for a scalar analogue of the pion 
form factor was studied with the help of the a-representation and the Mellin transformation. It was shown that, 
in the superrenormalizable (j)^ model, this diagram has the ln[Q 2 /m 2 ]/<3 4 behavior despite the fact that there 
is no logarithmic evolution in this model. The logarithm appears because the leading SD-pole 1/(J + 2) for the 
Mellin transform of this diagram can be obtained in two ways: from the small-pi integration (pl = cti + ot-i) and 
from the small-/9/{ integration (p^ = a 4 + as). There are no other possibilities. In particular, small-A integration 
(A = ot\ +a>2 +«3 + «4 + «5) gives a nonleading pole l/( J + 3). Hence, the leading term comes from a configuration 
in which the large momentum Q flows simultaneously through two subgraphs Vl = {01,172} and Vr = {04,05} 
while the momentum through the intermediate line 03 is small. Such a configuration was called in Ref. j^j the 
double-flow regime. 

In a renormalizable 3 6 ^ model the diagram shown in Fig. 11 has the \n 2 [Q 2 /m 2 ]/Q 4 behavior because the 
leading SD-pole 1/(J + 2) can be obtained in three ways: from small-A integration, from small- p^ integration and 
from small-pi? integration. The factorization for a scalar analog of the pion form factor in the model was 

studied in more detail in Ref. |6(|. It was shown there, in particular, that the In 2 Q 2 /m 2 behavior of the one-loop 
diagram results from the overlap of the evolution and the double-flow regime. In Ref. || , it was emphasized that 
the presence of the double-flow regime is a natural feature of exclusive ampitudes. Hence, to establish factorization, 
one should first check whether it is present or not. For the pion form factor in QCD (and other renormalizable 
models with spin-i quarks) its absence to all orders was demonstrated in Ref. ||. 

A rather peculiar double-flow contribution appears in a two-loop PQCD diagram for the nucleon form factors 
[^7). Its specifics is that it works for a term in which one takes only quark masses in the numerators of the 
propagators of the intermediate lines. Proceeding by a routine calculation, it is rather difficult to detect such a 
contribution among a wide variety of two loop terms. However, it is rather easy to find it if one has a guiding 
principle, such as the requirement that both the left and right components of a double-flow configuration should 
simultaneously give the leading power behavior. 



In our definitions of various distribution functions, we took the relevant matrix element and expressed it through 
an integral representation over the momentum fractions, incorporating the spectral condition < X < 1. Another 
approach is to introduce distribution functions by making a Fourier transform of the matrix element with respect 
to (pz) (cf. |^,^,^7)). One can easily derive the result of such a procedure by rewriting our representations in a 
form with a universal exponential in the r.h.s. Consider, e.g., the matrix element for the quark operator: 




■',p 2 )dr', 



IX. COMPARISON WITH OTHER APPROACHES AND NOTATIONS 



( P ',s'\$ a (0)zE(0,z;A)ip a (z)\p,s)\ z 



(9.1) 
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u(p',s')zu(p,s) I e -^ivz) F»(X;t)6(0<X <1)-Fl((-X;t)9(-1 + (<X <() 
-i+C L 



dX + "JC" 



The Fourier transformation would project out the function 

T£{X;t) = F£(X;t)6(0 <X< 1) - ^ c a (C - X; t) 6>( 1 + C < X < () (9.2) 

which a) coincides with the quark distribution for Q < X < \, b) reduces to the (minus) antiquark distribution for 
— 1 + C, < X < Q and c) is given by their difference for < X < (. The X- variable changes within the segment 
(— 1 + (, 1) centered at X — £/2, with the total range length equal to 2 — £. To avoid the nonsymmetric and 
(■-dependent limits, one can introduce the variable (cf. |10|) 

s _ X - C/2 
X ~ l-C/2 

which changes from —1 to 1. The ratio 



(9.3) 



(9.4) 



l-C/2 

is an alternative parameter characterizing the longitudinal momentum asymmetry of the nonforward matrix ele- 
ment. For t = and a massless hadron, it varies between and 2. The reversed relations are 

t _ j ~_ x + !j/2 ~ x ^/2 

C ~TT€72 ' TTIJ2 ■ X - C -T^j2- (9 - 5) 



Using translation invariance (cf. Eq. ( 4.15 )), one can easily derive that the operator with the quark fields taken 
at symmetric points —z/2, z/2 has a rather compact representation in terms of the i-variable: 

(p', s' | M-z/2)zE(-z/2, z/2; A)^ a {z/2) \p,s) |, 2=0 = u{p', s')zu( P , s) f e-"^H a (x, t) dx + «E a " , (9.6) 



where P = (p + p')/2 is the average momentum of the initial and final hadron (note, that (Pz) = (1 — C,/2)(pz) = 
(pz) I (l+£/2) ). This representation is equivalent to the definition of the off-forward parton distributions H a (x,^; f), 
E a (x,£,;t) introduced by X. Ji |l(]] (see also Jll|). Basically, the latter are related to our nonforward distributions 
by 

?£(X;t) = 0. + t/2)H a (x,Zit), (9.7) 

and similarly for other functions. The off-forward distributions H a (x,£,]t), etc. are defined both for positive and 
negative x. Depending on the value of x, one can distinguish three different components: quark (£/2 < x < 1), 
antiquark (— 1 < x < — £/2) and mixed "quark minus antiquark" (— ^/2 < x < ^/2) components of H. The mixed 
component corresponds evidently to the region < X < ( of the X- variable in which the nonforward distributions 
can be treated as distribution amplitudes. Since X(pz) = (x + £/2)(Pz) and (X — ()(pz) = (x — £/2)(Pz), the 
partons in this picture carry momenta (x + £/2)P and (x — £/2)P. Using Eqs.( |9.5| ),( |9.7|) , etc. one can convert our 
evolution kernels W^(X,Z) into those obtained earlier by X. Ji in Rcf. The gluonic matrix element can be 

also represented in the form of Eq. ([)!]) : 

(p'\z^Gl a (-z/2)E a \-z/2,z/2-A)G b au (z/2)\p)\ z , =f) 

= -u(p')zu(p){Pz) / e^ x( - Pz ^H g (i^;t)dx+^ . (9.8) 



2 



-l 



Due to the symmetry property H g (x,£;t) — H g (—x,£;t), integration over x in this case can be restricted to the 
< x < 1 region. Note, that in the forward limit £ = 0, t = 0, the function H g (x,£;t) reduces to xf g (x) (cf. Eq. 
(4.12)). To get an off-forward distribution reducing to f g (x), X. Ji [la] uses the definition equivalent to adding 
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a factor of x in the integrand on the right-hand side of Eq. (9.8): H g (x,£;t) — > xH^x^^t). However, x — 



corresponds to X — (/2 or to the middle-point Y — 1/2 of the distribution amplitude ^^(Y) (see Eq. ( |9.5| )), i.e., 
to a situation when the gluons carry equal fractions £p/2 of the original momentum p. Since H g (x,t;; t) is an even 
function of i, there are no evident reasons that it vanishes for x — 0. Hence, dividing H g (x, £; t) by x produces an 
artificial singularity of Hg l (x, £; t) for x = 0. 

Another parametrization for the nonforward matrix element of the gluon operator was proposed by Collins, 
Frankfurt and Strikman Jl~8| . Their definition of the nondiagonal gluon distribution f g (xi, x 2 ;t) is also based on 
the Fourier transformation. For positive values, their variables x\, X2 correspond to our fractions X and X — £ = X' , 
respectively. In our notations, the function f g {x\ — X, x 2 = X — £; t) can be written as 

f g (X, X-C,t)= ?1{X; t) . (9.9) 

The factor 1/X(X — £) was motivated by the necessity to cancel the inverse factor which may emerge from the 
derivatives present in the field-strength tensor G^ v . Actually, this expectation is not supported by perturbative 
calculations. Take, e.g., the evolution kernel P 9 ®(X,Z). It can be treated as a perturbative, leading logarithm 



approximation for the gluon distribution inside a quark (cf. [HJ). According to Eq. ( 5.28 ), Pf^{X,Z) does not 



vanish for X = £. If T 9 (X) does not vanish for X = (, the function T 9 {X; t) does not vanish both for X = and 
X = £ and f g (xi, x 2 ) is singular both for x\ = and x 2 = 0. 

In fact, the combination T 9 {X)/{X — ie)(X — £ + ie) appears in our parametrization ( [7. 14 ) for the matrix 



element of the operator constructed from two vector potentials A v taken in the light-cone gauge. In this sense, 
f g {xi, x 2 ) or, what is the same, J- 9 (X)/X(X — £) can be treated as a basic gluon distribution given by the matrix 
element of the product of fundamental gluonic fields A^A V rather than by that of the secondary fields G^pG?. 
Note, however, that if f g (xi,x 2 ), i.e. T 9 (X)/X(X ~ £), has no singularities, then the meson electroproduction 
amplitude has no imaginary part at leading twist. Since this is impossible, f g (xi,x 2 ) must have singularities, and 
one may wish to explicitly display them specifying their nature, e.g., 1/(^2 — ie), l/(^i + *e)- This goal is achieved 
automatically if F®{X) is used as the basic distribution. 

In our approach, the starting point is the double distribution F g (x,y;t) defined through the nonforward matrix 
element of the gauge-invariant gluonic operator 

(p'\z^G a , a (0)E ab (0,z;A)G h a ^)\p)\^o 

= u(p')zu(p)(pz) dx i (z-ix(pz)-iy(rz) + e ix(pz)-iy(rz)^ F g {x , y; t) 9(x + y < I) dy . (9.10) 

As explained earlier, in perturbation theory the spectral properties < {x, y,x + y} < 1 can be proved to any order 
with the help of the a- representation. Furthermore, the function F g (x,y;t) does not depend on the ^-parameter. 
The family of £-dependent nonforward gluon distributions T®(X; t) is obtained from F g {x, y; t) by integration over 



y (see (2.46)): 



,min{X/C,X/C} 

T 9 AX;t)= F g (X-y(,y;t)dy. (9.11) 

Jo 

Recall that the double distribution F g {x,y;t) can be treated as a distribution function with respect to x and 
as a distribution amplitude with respect to y. This physical interpretation suggests that F g {x,y;t) is a regular 
function for all values of y and for at least nonzero values of x. We made this reservation because the evolution 
asymptotically makes F g (x,y;t; /i) (we added the dependence on the factorization scale /i) proportional to 8(x) 
as /i — > oo. In this situation, F g (x,y;t; fi) is singular at x = 0. However, the ^(x)-term still produces a regular 
nonforward distribution T 9 {X;t), though confined to the restricted region < X < (\ 

Assuming that the double distribution F g (x,y;t; /i) is finite everywhere, we conclude that the nonforward dis- 
tribution F 9 (X;t;fi) in this case is also finite for all < X < 1 and, moreover, that it vanishes for X = 0. As 
discussed earlier, the latter property is vital for factorization. If it is not fulfilled, the A-integral in the lowest-order 
expression diverges at the end-point X — 0, where the 1/(A — ie) prescription is of no help. One may think 
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that this problem can be avoided if one uses the function J-^{X\t) defined through the Fourier transformation 
with the variable X changing from — 1 + £ to 1. Since the point X = is inside the integration region, the 
1/(X — ie) prescription apparently may help. Note, however, that if our function T^{X]t) does not vanish for 
X = 0, the Fourier transform T^(X; t) is not continuous both for X — and X — (. As a result, the singularities 
of !Ff(X; t)/(X — ie)(X — C + ie) are not integrable. 



X. CONCLUSIONS 



In this paper, we discussed basic properties of nonforward parton distributions, a new type of functions accu- 
mulating nonperturbative information about hadron dynamics. We demonstrated that there are two basic ways to 
describe asymmetric matrix elements (p'\O(0, z)\p) of quark and gluon light-cone operators O(0, z). One possibil- 
ity is to introduce double distributions F(x, y; i) which are independent of the longitudinal momentum asymmetry 
(or skewedness) parameter £ — 1 — (p'z)/(pz) of the matrix element and refer to the light-cone fractions xp, yr 
of the original hadron momentum p and momentum transfer r = p' — p carried by the active parton. Another 
approach is to use nonforward distribution functions J-^(X;t) which specify the light-cone projection of the total 
momentum Xp = xp + yr carried by the parton. These functions T(^{X]i) explicitly depend on Both types of 
distributions have hybrid properties, in some aspects resembling usual parton distribution functions and in other 
ones the distribution amplitudes. Their ^-dependence is analogous to that of hadronic form factors. The use of 
!F^(X;t) is more convenient for ultimate applications to hard PQCD processes, resulting in a formalism that is 
very similar to the standard PQCD parton picture. On the other hand, the double distributions F(x,y;t) have 
more transparent spectral properties which has serious advantages at the foundation stages of the PQCD analy- 
sis. In this paper, we concentrated on general aspects of the theory of nonforward distributions and their uses. 
There are many interesting applications to deeply virtual Compton scattering and hard exclusive electroproduction 
processes which require further, more specific studies of the nonforward distribution functions including modeling 
their nonperturbative low-energy shape, logarithmic PQCD evolution, calculation of nonlogarithmic higher-order 
corrections, etc. Work in this direction has already been started [^0|- p3[ . 
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FIG. 1. Scalar model analogs of a),b) virtual forward Compton amplitude and c),d) deeply virtual Compton scattering. 
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FIG. 2. Longitudinal momentum flow for two components of the asymmetric distribution function T^X): a) X > £ and 

b)x<c- 
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FIG. 3. a) Structure of momentum integral defining the asymmetric distribution function (X) . b) Cut of parton-hadron 
amplitude corresponding to the residue for the region X > (. c) Cut of parton-hadron amplitude corresponding to the residue 
for the region X < 
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FIG. 8. a) General structure of the leading SD contribution to the DVCS amplitude in QCD. b) SD configuration with 
two long-distance parts. c) Matrix element with double-logarithmic dependence on the UV cut-off parameter fi. 




FIG. 9. Hard exclusive meson electoproduction process: a) Leading SD-contribution with quark nonforward distribution; 
b) Soft contribution. 
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FIG. 10. Gluon contribution to hard exclusive meson electoproduction amplitude. 
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FIG. 11. Double-flow regime for the scalar analog of a meson form factor. 
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